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PsJ Abstract The MEG (Mu to Electron Gamma) experiment Keywords High resolution time measurement; Rare muon 
has been running at the Paul Scherrer Institut (PSI), Switzer- decays; Magnetic spectrometer; Liquid Xenon calorimeter; 
land since 2008 to search for the decay yU^ e^y by us- Lepton Flavour Violation 
ing one of the most intense continuous /u^ beams in the 
I world. This paper presents the MEG components: the posi- 
J> Iron spectrometer, including a thin target, a superconducting 

magnet, a set of drift chambers for measuring the muon de- Contents 
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Fig. 1 Schematic view of the MEG detector showing one simulated signal event emitted from the target. 



1 Introduction 

A search for the Charged Lepton Flavour Violating (CLFV) 
decay — > e^y, the MEG experiment (see [ IJ and refer- 
ences therein for a detailed report of the experiment moti- 
vation, design criteria and goals) is in progress at the Paul 
Scherrer Institut (PSI) in Switzerland. Preliminary results 
have already been published |I2]|3]- The goal is to push the 
sensitivity to this decay down to ~ 5x10"'-' improving the 
previous limit set by the MEGA experiment, 1 .2 x 10"" L4J, 
by a factor 20. 

CLFV processes are practically forbidden in the Stan- 
dard Model (SM), which, even in presence of neutrino masses 
and mixing, predicts tiny branching ratios (BR <sc 10"^") 
for CLFV decays. Detecting such decays would be a clear 
indication of new physics beyond the SM, as predicted by 
many extensions such as supersymmetry [ 5_ | . Hence, CLFV 
searches with improved sensitivity either reveal new physics 
or constrain the allowed parameter space of SM extensions. 

In MEG positive muons stop and decay in a thin target 
located at the centre of a magnetic spectrometer The signal 
has the simple kinematics of a two-body decay from a parti- 
cle at rest: one monochromatic positron and one monochro- 
matic photon moving in opposite directions each with an 
energy of 52.83 MeV (half of the muon mass) and being 
coincident in time. 

This signature needs to be extracted from a background 
induced by Michel (/i^ — > e^vv) and radiative (//^ — > e^yvv) 
muon decays. The background is dominated by accidental 
coincidence events where a positron and a photon from dif- 
ferent muon decays with energies close to the kinematic 



limit overlap within the direction and time resolution of the 
detector. Because the rate of accidental coincidence events 
is proportional to the square of the yu^ decay rate, while 
the rate of signal events is proportional only to the yu^ de- 
cay rate, direct-current beams allow a better signal to back- 
ground ratio to be achieved than for pulsed beams. Hence 
we use the PSI continuous surface muon beam with inten- 
sity ~ 3 X lOV^/s (see Sect.|2ll. 

A schematic of the MEG apparatus is shown in Fig. [1] 
A magnet, COBRA (COnstant Bending RAdius), generates 
a gradient magnetic field, for the first time among particle 
physics experiments, with the field strength gradually de- 
creasing at increasing distance along the magnet axis from 
the centre. 

This configuration is optimised to sweep low-momentum 
positrons from Michel decays rapidly out of the magnet, and 
to keep the bending radius of the positron trajectories only 
weakly dependent on their emission angle within the accep- 
tance region (see Sect.O. 

The positron track parameters are measured by a set of 
very low mass Drift CHambers (DCH) designed to minimise 
the multiple scattering (see Sect. |4|i. The positron time is 
measured by a Timing Counter (TC) consisting of scintil- 
lator bars read out by PhotoMultiplier Tubes (PMT) (see 
Sect.|5j. 

For y-ray detection, we have developed an innovative 
detector using Liquid Xenon (LXe) as a scintillation mate- 
rial viewed by PMTs submersed in the liquid. This detector 
provides accurate measurements of the y-ray energy and of 
the time and position of the interaction point (see Sect.|6]l. 



3 



We monitor the performance of the experiment contin- 
uously by introducing a variety of calibration methods (see 
Sect. IT). 

We employ a flexible and sophisticated electronics sys- 
tem for triggering the fi^ — > e^y candidate events along with 
calibration data (see Sect.[8]l. 

Special considerations influenced the development of the 
data acquisition system. We record signals from all detec- 
tors with a custom-designed waveform digitiser, the Domino 
Ring Sampler (DRS) (see Sects.l9]and[T0t. 

We evaluate detector response functions using data with 
help of detector simulations (see Sect.fTTTi. 

In this paper, we mostly refer to a cylindrical coordinate 
system (r, 0, z) with origin in the centre of COBRA. The z- 
axis is parallel to the COBRA axis and directed along the in- 
coming muon beam. The axis defining <p - (the x-axis of 
the corresponding Cartesian coordinate system) is directed 
opposite to the centre of the LXe detector (as a consequence 
the y-axis is directed upwards). Positrons move along trajec- 
tories of decreasing (p coordinate. When required, the polar 
angle with respect to the z-axis is also used. The region 
with z < is referred as upstream, that with z > as down- 
stream. 

The angular acceptance of the experiment is defined by 
the size of the LXe fiducial volume, corresponding approxi- 
mately to e (jTT, ^n) and | cos 6\ < 0.35, for a total accep- 
tance of ~ 11%. The spectrometer is designed to fully accept 
positrons from jj^ — > e^y decays when the y-ray falls in the 
acceptance region. 



2 Beam Line 

2.1 Experimental Requirements 

In order to meet the stringent background requirements nec- 
essary for a high-sensitivity, high-rate coincidence experi- 
ment such as MEG, one of the world's highest intensity con- 
tinuous muon beams, the 7rE5 channel at PSI, is used in con- 
junction with the high-performance MEG beam transport 
system which couples the nE5 channel to the COBRA spec- 
trometer (see Sect. 13). This allows more than 10** /^Vs, with 
an optimal suppression of beam-correlated background, to 
be transported to the ultra-thin stopping-target at the centre 
of COBRA. The main requirements of such a low-energy 
surface muon beam are: high-intensity with high-transmission 
optics at 28 MeV/c, close to the kinematic edge of stopped 
muon decay and to the maximum intensity of the beam mo- 
mentum spectrum; small beam emittance to minimise the 
stopping-target size; momentum-biteQof less than 10% with 
an achromatic final focus to provide an almost monochro- 

' Selected range of momenta transmitted by the channel. 



matic beam with a high stopping density |j and enabling 
muons to be stopped in a minimally thick target; minimi- 
sation and separation of beam-related background, such as 
beam originating from ;r"-decay in the production target 
or decay particles produced along the beam line; minimisa- 
tion of material interactions along the beam line, requiring 
both vacuum and helium environments to reduce the mul- 
tiple scattering and the photon-production probability for 
annihilation-in- flightd or Bremsstrahlung. 

Furthermore, the experiment requires an arsenal of cali- 
bration and monitoring tools, including, on the beam side, 
the possibility of sources of monochromatic photons and 
positrons for calibration purposes. This in turn requires a 
70.5 MeV/c stopped negative-pion beam tune to produce 
photons from the charge-exchange {n^p — > Tv^n) and ra- 
diative capture reactions (n^p — » y«) on a liquid hydro- 
gen target (see Sect. 17.3.2b . as well as a 53 MeV/c positron 
beam tune to produce monochromatic Mott-scattered posi- 
trons (see Sect. 17.2b . 

Finally, the MEG beam transport system also allows the 
introduction of external targets and detectors, as well as a 
proton beam line, from a dedicated Cockcroft- Walton (C- 
W) accelerator (see Sect. 17.3. IL to the centre of the spec- 
trometer helium volume. This is achieved without affecting 
the gas environment, by means of specially designed end- 
caps. The protons from C-W accelerator interact with a ded- 
icated target to produce y-ray lines induced by nuclear reac- 
tions. 

2.2 7rE5 Beam line and MEG Beam Transport System 

The 7rE5 channel is a 165° backwards-oriented, windowless, 
high-acceptance (150 msr), low-momentum (< 120 MeV/c), 
dual-port tt , yu" or e-channel. For the MEG experiment, 
the channel is tuned to H-28 MeV/c with a momentum-bite 
of between 5-7% FWHM, depending on the opening of 
the momentum selecting slits placed in the front-part of the 
channel. This constitutes an optimal high-intensity surface 
muon tune (stopped pion-decay at the surface of the pro- 
duction target) at backward pion production angles [(6]|7]. 
A schematic of the 7rE5 channel and the connecting MEG 
beam line is shown in Fig. |2] The main production target 
is a rotating radiation-cooled graphite wheel of 4 cm thick- 
ness in the beam direction. It is connected to the 7rE5 area 
by a quadrupole and sextupole channel that couples to the 
initial extraction and focussing elements of the MEG beam 
line, quadrupole triplet I, which exits the shielding wall and 
provides optimal high transmission through the Wien-filter 

measure of the muon stopping rate per unit target volume, influ- 
enced by both the beam momentum and bite as well as the target ma- 
terial properties. 

'The reaction e^e" yy with e* from /j* Michel decay and e" in the 
material. 
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Fig. 2 Schematic of the ;rE5 Channel and the MEG beam transport system showing the beam elements and detectors: Wien-filter, collimator system 
and the superconducting Beam Transport Solenoid (BTS) and central degrader system. Coupled to these are the He-filled COBRA spectrometer, 
consisting of a gradient-field magnet and central target system. The downstream insertion system allowing the remote insertion of the Cockcroft- 
Walton beam pipe, for calibration purposes, is also shown. 



(E X B field separator). With an effective length of about 
83 cm for both the electric and magnetic fields, together 
with a 200 kV potential across the 19 cm gap of the elec- 
trodes, a mass separation equivalent to an angular separation 
of +88 mrad and -25 mrad between muons and positrons 
and between muons and pions is achieved, respectively. This 
together with Triplet II and a collimator system (equivalent 
to 1 1 Xo of Pb), placed after the second triplet, separates 
the muons from the eight-times higher beam positron con- 
tamination coming from the production target. A separation 
quality between muons and beam positrons of 8.1 cr^ is so 
achieved, corresponding to a 12 cm physical separation at 
the collimator system as shown in Fig. [3] This allows an al- 
most pure muon beam0to propagate to the Beam Transport 
Solenoid (BTS) (see Sect. 12.31 ), which couples the injection 
to the superconducting COBRA spectrometer in an iron-free 
way. By means of a 300 fim thick Mylar ® degrader system 
placed at the intermediate focus, the BTS also minimises the 
multiple scattering contribution to the beam and adjusts the 
muon range for a maximum stop-density at the centre of the 
muon target, placed inside a helium atmosphere, at the cen- 
tre of the COBRA spectrometer. 

2.3 The Beam Transport Solenoid (BTS) 

The BTS, a 2.8 m long iron-free superconducting solenoid, 
with a 38 cm warm-bore aperture, coupling directly to the 
beam line vacuum, is shown in Fig.|4] With a maximal cur- 
rent of 300 A a field of 0.54 T can be reached, higher than 
the 0.36 T required for an intermediate focus at the central 
degrader system for muons. The double-layered 2.63 m long 

''Positron suppression factor ~ 850 under the muon peak. 



coil is made up of a 45% Ni/Ti, 55% Cu matrix conductor 
with an inductance of 1.0 H so yielding a maximum stored 
energy of 45 kJ. 

On the cryogenic side, a dedicated Uquid helium (LHe) 
port to the central cryoplant is used to transport the super- 
critical helium at ~ 10 bar to two Joule-Thompson valves 
(in the upper tower in Fig. |4]i, which are used to liquefy and 
automatically regulate the level of the 150 of LHe in the 
cryostat. A separate liquid nitrogen line cools an outer heat- 
shield used to minimise the LHe losses. 

Furthermore, three iron-free correction dipole magnets 
(one cos 0-type placed upstream of the BTS) one horizon- 
tal and one vertical conventional-dipole type, both mounted 
on the downstream side of the BTS, are employed to ensure 
an axial injection into both solenoids. The three dipole mag- 
nets compensate for a radial asymmetry in the fringe fields 
induced by the interaction of the stray field from the large- 
aperture, iron-free COBRA magnet with an iron component 
of the hall-floor foundations. An example of the optics in 
Fig. |5] shows the 2cr beam envelopes both vertically and 
horizontally along the 22 m long beam line, together with 
the apertures of the various elements and the dispersion tra- 
jectory of particles having a 1% higher momentum (dotted 
line). 



2.4 Spectrometer Coupling System 

The COBRA spectrometer is sealed at both ends by end- 
caps which constitute the final part of the beam line, sep- 
arating the vacuum section from the helium environment 
inside. These specially designed 1160 mm diameter light- 
weight (35^0 kg), only 4 mm-thick aluminium constructs 
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Fig. 3 Separator scan plot, measured post collimator system. The black 
dots represent on-axis, low discriminator threshold intensity measure- 
ments at a given separator magnetic field value during the scan, for a 
fixed electric field value of -195 kV. The red curve represents a double 
Gaussian fit to the data points, with a constant background. A separa- 
tion of 8.1 muon beam cr^ is found, corresponding to more than 12 cm 
separation of beam-spot centres at the collimator The raw muon peak 
measured at low discriminator threshold also contains a contribution 
from Michel positrons due to decays in the detector These can eas- 
ily be distinguished by also measuring at high threshold, where only 
muons are visible on account of their higher energy loss and hence 
higher pulse-height. 



minimise the material available for possible photon back- 
ground production, crucial to a background-sensitive exper- 
iment. They also integrate the beam vacuum window on the 
upstream side and an automated thin-windowed (20 fim) 
bellows insertion system, capable of an extension-length of 
about 1.7 m. The bellows insertion system is used to intro- 
duce various types of targets such as a liquid hydrogen target 
used in conjunction with the pion beam, or a lithium tetrab- 
orate photon production target used together with the C-W 
accelerator, on the downstream side. 



2.5 Commissioning 

The commissioning of the MEG beam line was completed 
with the first "engineering run" in late 2007. During com- 
missioning, several beam optical solutions were studied, in- 
cluding the investigation of both nE5 branches ("U"-branch 
and "Z"-branch) as well as different beam-coiTelated back- 
ground separation methods: a combined degrader energy- 
loss/magnetic-field separation method or a crossed (E x B)- 
field separator (Wien filter). The optimum solution favoured 
the "Z"-channel with a Wien filter, both from the experimen- 




Fig. 4 The Beam Transport Solenoid (BTS), shown coupled to the 
beam line between Triplet II and the COBRA spectrometer The main 
cryogenic tower with the current leads and the loule-Thompson valve 
tower with the LHe transfer line from the cryoplant coupled to it as well 
as the smaller LN2 transfer line in the background can be seen. The two 
coiTection dipole magnets on the downstream side are not shown here. 



tal layout and the rate and background-suppression point of 
view. 

In order to maximise the beam transmission rate and 
the stopping density in the target, an understanding of the 
momentum spectrum and channel momentum acceptance 
is necessary. This was carefully studied by measuring the 
muon momentum spectrum in the vicinity of the kinematic 
edge of stopped pion decay (29.79 MeV/c). As with all phase- 
space measurements performed at high rate, a scanning tech- 
nique was used: either one using a small cylindrical "pill" 
scintillator of 2 mm diameter and 2 mm thickness coupled 
to a miniature PMT |8 1, or one using a thick depletion layer 
Avalanche PhotoDiode (APD), without a scintillator, in the 
case of measurements in a magnetic field. Both these so- 
lutions gave an optimum pulse-height separation between 
positrons and muons. This allows integrated beam-spot in- 
tensity measurements up to 2 GHz to be made. 

For the momentum spectrum, taken with the full 7% 
FWHM momentum-bite and shown in Fig.[6l the whole beam 
line was tuned to successive momentum values and the inte- 
gral muon rate was then determined by scanning the beam- 
spot using the beam line as a spectrometer This allows to 
determine both the central momentum and the momentum- 
bite by fitting a theoretical distribution, folded with a Gaus- 
sian resolution function corresponding to the momentum- 
bite, plus a constant cloud-muon background over this Um- 
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Fig. 5 2" -order TRANSPORT calculation showing the vertical (top-half) and horizontal (bottom-half) beam envelopes plotted along the beam- 
axis from close to target E up to the centre of COBRA, a distance of 22 m. The vertical and horizontal scales differ 33:1. The coloured elements 
show the vacuum chamber apertures of the various beam elements, while the dotted line represents the central dispersion trajectory for particles 
with a 1 % higher momentum. 



ited momentum range. The red curve corresponding to the 
fitted function shows the expected empirical behaviour 
below the kinematic edge. The cloud-muon content below 
the surface-muon peak, which has its origin in pion decay- 
in-flight in the vicinity of the production target, was deter- 
mined to be 3.3% of the surface-muon rate. This was de- 
termined from negative muon measurements at 28 MeV/c 
and using the known pion production cross-section differ- 
ences ||9], since negative surface-muons do not exist due to 
the formation of pionic atoms and hence the small cloud- 
muon contribution is much easier to measure in the absence 
of a dominating surface-muon signal. 

Another important feature investigated was the muon po- 
larisation on decaying: since the muons are born in a 100% 
polarised state, with their spin vector opposing their mo- 
mentum vector, the resultant angular distribution on decay 
is an important issue in understanding the angular distribu- 
tion of the Michel positrons and the influence of detector 
acceptance. The depolarisation along the beam line, influ- 
enced by divergence, materials, magnetic fields and beam 
contaminants was studied and incorporated into the MEG 
Monte Carlo, which was used to compare results with the 
measured angular distribution asymmetry (upstream/down- 
stream events) of Michel positrons. The predicted muon po- 
larisation of 91% was confirmed by these measurements, 
which yielded a value of 89 ± 4 % HlOlll II . The main single 
depolarising contribution is due to the mean polarization of 
the 3.3% cloud-muon contamination, which opposes that of 
the surface muons. However, because the depolarisation in 
the polyethylene target is totally quenched by the strength of 
the longitudinal magnetic field in COBRA, a high polariza- 
tion can be maintained [.lOJ . 
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Fig. 6 Measured MEG muon momentum spectrum with fully open 
slits. Each point is obtained by optimising the whole beam line for the 
coiTesponding central momentum and measuring the full beam-spot 
intensity. The red line is a fit to the data, based on a theoretical P^-^ 
behaviour, folded with a Gaussian resolution function corresponding 
to the momentum-bite plus a constant cloud-muon background. 



In order to achieve the smallest beam-spot size at the 
target, the position of the final degrader as well as the polar- 
ity combination of COBRA and the BTS were studied. The 
smallest beam divergences and hence the smallest beam- 
spot size at the target were achieved by having a reverse 
polarity magnetic field of the BTS with respect to COBRA. 
This is due to the fact that the radial field is enhanced at 
the intermediate focus between the two magnets, thus effec- 
tively increasing the radial focussing power and minimis- 
ing the growth of the beam divergences through the CO- 
BRA volume. The final round beam-spot at the target, with 
o".i,>' ~ 10 mm is shown in Fig.|7l 
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Fig. 7 Scanned horizontal and vertical muon beam-spot profiles mea- 
sured at the centre of COBRA. The red lines show the Gaussian fits 
to the distributions that yield cr,,. ~ 10 mm and a total stopping in- 
tensity of ~ 3 X 10''ju+/s at a proton beam intensity of 2.2mA and 
momentum-bite of 5.0% FWHM. 



2.6 Beam line Performance 

Although the beam line was tuned for maximal rate with 
minimal positron beam contamination, the final beam stop- 
ping intensity was determined by a full analysis of the de- 
tector performance and an evaluation of its sensitivity to the 
//^ — > e^y decay. This yielded an optimal muon stopping 
rate in the target of 3 x 10^ //^/s. 

Table 1 MEG Final Beam Line Characteristics 



Parameter 


Value 


Solid Angle zifl 


150 msr 


Beam Line length L 


22 m 


Muon Central Momentum 


28 MeV/c 


Pion Central Momentum P;, 


-70.5 MeV/c 


Positron Central Momentum P, 


53 MeV/c 


Muon Momentum-bite AP/P 


5.0% (full-width) 


Beam-spot size ctx.y 


9-10 mm 


beam e^/ju^ ratio (Sep. off) 


8 


ju*/e* separation (Sep. on) 


8.1o"p, 12 cm 


Optimal )X* stopping intensity 


3x lOV/s 



Furthermore, for both the pion beam and Mott positron 
beam calibration tunes, mentioned in Sect. 12.11 which in- 
volve prompt particle production at the proton target, the 
central momentum was chosen to give a maximal time-of- 



flight separation between the wanted particles and the un- 
wanted background particles at the centre of COBRA, based 
on the use of the accelerator modulo 19.75 ns radio-frequency 
signal as a clock, which has a fixed time relationship to the 
production at the proton target. In the case of the Mott po- 
sitron beam the final momentum should also be close to the 
signal energy of 52.83 MeV/c. A summary of the final beam 
parameters is given in Table [1] 



2.7 Target 

2.7.7 Requirements 

Positive muons are stopped in a thin target at the centre of 
COBRA. The target is designed to satisfy the following cri- 
teria: 

1 Stop > 80 % of the muons in a limited axial region at the 
magnet centre. 

2 Minimise conversions of y-rays from radiative muon de- 
cay in the target. 

3 Minimise positron annihilation-in-flight with y-rays in 
the acceptance region. 

4 Minimise multiple scattering and Bremsstrahlung of po- 
sitrons from muon decay. 

5 Allow the determination of the decay vertex position and 
the positron direction at the vertex by projecting the re- 
constructed positron trajectory to the target plane. 

6 Be dimensionally stable and be remotely removable to 
allow frequent insertion of ancillary calibration targets. 

Criterion 1 requires a path-length along the incident muon 
axial trajectory of ~ 0.054 g/cm^ for low-Z materials and for 
the average beam momentum of 19.8 MeV/c (corresponding 
to a kinetic energy of 1.85 MeV) at the centre of COBRA. 
To reduce scattering and conversion of outgoing particles, 
the target is constructed in the form of a thin sheet in a ver- 
tical plane, with the direction perpendicular to the target at 
6 ~ 70°, oriented so that the muon beam is incident on the 
target side facing the LXe detector (see Fig.[T]i. The require- 
ment that we be able to deduce the positron angle at the ver- 
tex sets the required precision on the knowledge of the posi- 
tion of the target in the direction perpendicular to its plane. 
This derives from the effect of the track curvature on the in- 
ferred angle; an error in the projected track distance results 
in an error in the positron angle at the target due to the curva- 
ture over the distance between the assumed and actual target 
plane position. This angle error is larger when the positrons 
are incident at a large angle with respect to the normal vec- 
tor; for an incidence angle of 45°, an error in the target plane 
position of 500 fim results in an angle error of ~ 5 mrad. 
This sets the requirements that the target thickness be below 
500 fj.m (eliminating low-density foam as a target material) 
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and that its position with respect to the DCH spectrometer 
be measured to ~ lOOjt/m. 

2.7.2 Realisation 

To meet the above criteria, the target was implemented as a 
205 jjm thick sheet of a low-density (0.895 g/crn^) layered- 
structure of a polyethylene and polyester with an elliptical 
shape and semi-major and semi-minor axes of 10 cm and 
4 cm. The target was captured between two rings of Roha- 
cell ® foam of density 0.053 g/cm^ and hung by a Delrin ® 
and Rohacell stem from a remotely actuated positioning de- 
vice attached to the DCH support frame. To allow a software 
verification of the position of the target plane, as described 
in Sect. 12.7.51 eight holes of radius 0.5 cm were punched in 
the foil. 

2.7.3 Automated insertion 

The positioning device consists of a pneumatically oper- 
ated, bi-directional cylinder that inserts and retracts the tar- 
get. Physical constraints on the position of the cylinder re- 
quired the extraction direction to not be parallel to the target 
plane; hence the position of the plane when inserted was 
sensitive to the distance by which the piston travelled. This 
distance was controlled by both a mechanical stop on the 
piston travel and a sensor that indicated when the target was 
fully inserted. 

2. 7.4 Optical alignment 

The position of the target plane is measured every run period 
with an optical survey that images crosses inscribed on the 
target plane; this survey is done simultaneously with that 
of the DCH. The relative accuracy of the target and DCH 
positions was estimated to be ±(0.5, 0.5, 0.8) mm in the (x, 
y, z) directions. 

2.7.5 Software Alignment 

The optical survey was checked by imaging the holes in the 
target. We define z' as the position parallel to the target plane 
and perpendicular to the y coordinate and x' as the position 
perpendicular to the target plane. The y and z' positions are 
easily checked by projecting particle trajectories to the nom- 
inal target plane and determining the z' and y hole coordi- 
nates by fitting the illuminations to holes superimposed on 
a beam profile of approximately Gaussian shape. We note 
that these target coordinates are not important to the anal- 
ysis, since a translation of the target in the target plane has 
no effect on the inferred vertex position. The x' coordinate 
was checked by imaging the y positions of the holes as a 
function of the track direction at the target (p^. Figure |8] il- 
lustrates the technique. The x' target position is determined 



by requiring that the target hole images be insensitive to (pg. 
Figure |9] shows the illumination plot for the 20 11 data sam- 
ple. Figure[TO]shows the normalised y distribution of events 
with 20° < (pe < 30° in a shce at fixed z' centred on a target 
hole (in red) and in adjacent slices (in blue). The difference 
(in black) is fitted with a Gaussian. Figure [TT] shows the ap- 
parent y position of one of the target holes versus (pe, fitted 
with a tan(0(,) function, from which we infer that the x' offset 
with respect to the optical survey is 0.2+0.1 mm. By imag- 
ing all the holes (i.e. finding x'{y,z')) both the planarity and 
the orientation of the target plane was determined. 

Y A assumed target 




► 

X' 

Fig. 8 An illustration of the target hole imaging technique. If there is 
an offset in the x' direction between the assumed target position and the 
true target position, the apparent v position of a target hole will depend 
on the angle of emission at the target (p^. The shift in y is approximately 
proportional to the shift in the x' direction by ianiipe)- 




Fig. 9 The target illumination plot for the 20 1 1 data sample. 



2.7.6 Multiple scattering contribution from the target 

The target itself constitutes the most relevant irreducible con- 
tribution to the positron multiple scattering and therefore to 
the angular resolution of the spectrometer. 
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Fig. 10 The nonnalised y distribution of events with 20° < ipe < 30° 
in a slice at fixed z' centred on a target hole (in red) and in adjacent 
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Fig. 11 The apparent y position of a target hole versus (p^ fitted with a 
tan(0(,) function, from which we infer that the x' off'set with respect to 
the optical survey is 0.2±0.1 mm. 
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Fig. 12 Distribution of the /i+ decay vertex along x' . 



The simulated distribution of the ju^ decay vertex along 
x' is shown in Fig. [12] The average is slightly beyond the 
centre of the target. Averaging over the decay vertex posi- 
tion and over the positron direction in the MEG acceptance 
region, the difference between the positron direction at the 



yU^ decay vertex and that at the target boundary follows a 
Gaussian distribution with cr - 5.2 mrad. 



3 COBRA magnet 

3.1 Concept 

The COBRA (COnstant Bending RAdius) magnet is a thin- 
wall superconducting magnet generating a gradient mag- 
netic field which allows stable operation of the positron spec- 
trometer in a high-rate environment. The gradient magnetic 
field is specially designed so that the positron emitted from 
the target follows a trajectory with an almost constant pro- 
jected bending radius weakly dependent on the emission po- 
lar angle 9 (Fig.[T3]'a)). Only high-momentum positrons can 
therefore reach the DCH placed at the outer radius of the 
inner bore of COBRA. Another good feature of the gradi- 
ent field is that the positrons emitted at cos 9 ~ are quickly 
swept away (Fig. [T3l' b)). This can be contrasted to a uniform 
solenoidal field where positrons emitted transversely would 
turn repeatedly in the spectrometer The hit rates of Michel 
positrons expected for the gradient and uniform fields are 
compared in Fig. [14] and indicate a significant reduction for 
the gradient field. 

The central part of the superconducting magnet is as thin 
as 0.197 X() so that only a fraction of the y-rays from the 
target interacts before reaching the LXe detector placed out- 
side COBRA. The COBRA magnet (see Fig.[T5]) is equipped 
with a pair of compensation coils to reduce the stray field 
around the LXe detector for the operation of the PMTs. 

The parameters of COBRA are listed in Table|2] 



3.2 Design 

3.2.1 Superconducting magnet 

The gradient magnetic field ranging from 1.27 T at the cen- 
tre to 0.49 T at either end of the magnet cryostat is gen- 
erated by a step structure of five coils with three diff'erent 
radii (one central coil, two gradient coils and two end coils) 
(Fig. [T6l l. The field gradient is optimised by adjusting the 
radii of the coils and the winding densities of the conductor 
The coil structure is conductively cooled by two mechan- 
ical refrigerators attached to the end coils; each is a two- 
stage Gifford-McMahon (GM) refrigerator [12] with a cool- 
ing power of 1 W at 4.2 K. The thin support structure of the 
coil is carefully designed in terms of the mechanical strength 
and the thermal conductivity. The basic idea is to make the 
coil structure as thin as possible by using a high-strength 
conductor on a thin aluminium mechanical support cylinder. 
Figure[T7]shows the layer structure of the central coil, which 
has the highest current density. A high-strength conductor 
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Table 2 Parameters of the COBRA magnet. 

Coil 



Central Gradient Inner end Outer end Compensation 



Conductivity 


Super 


Super 


Super 


Super 


Resistive 


Inner dia. (mm) 


699.1 


809.1 


919.1 


919.1 


2210 


Outer dia. (mm) 


711.6 


820.6 


929.5 


929.5 


2590 


Length (mm) 


240.3 


110.4 


189.9 


749.2 


265 


z-coordinate of coil centre(mm) 


0.8 


±235 


±405.4 


±874.95 


±1190 


Layers 


4 


4 


3 


3 


14 


Turns (total) 


1068 


399 


240 


1548 


280 


Winding density(Tums/m) 


4444.4 


3614.1 


1263.8 


2066.2 


1056.6 


Inductance(H) 


1.64 


0.62 


0.35 


2.29 


0.54 


Current (A) 


359.1 


359.1 


359.1 


359.1 


319.2 


Energy E (kj) 


106 


40 


23 


148 


35 


Weight M (kg) 


9 


4 


7 


28 


1620 


EIM (kJ/kg) 


11.8 


10.0 


3.3 


5.3 


0.02 




(a) 




(b) 



Fig. 13 Concept of the gradient magnetic field of COBRA. The posi- 
trons follow trajectories at a constant bending radius weakly dependent 
on the emission angle 9 (a) and those transversely emitted from the tar- 
get (cos 6 ~ 0) are quickly swept away from the DCH (b). 



is wound in four layers inside the 2mm-thick aluminium 
support cylinder Pure aluminium strips with a thickness of 
100 //m are attached on the inner surface of the coil structure 
in order to increase the thermal conductivity. Several quench 
protection heaters are attached to all coils in order to avoid 
a local energy dump in case of a quench. The high thermal 
conductivity of the coil structure and the uniform quench in- 
duced by the protection heaters are important to protect the 
magnet. 
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Fig. 14 Hit rate of the Michel positrons as a function of the radial 
distance from the target in both the gradient and uniform field cases. 
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Fig. 15 Cross-sectional view of the COBRA magnet. 

A high-strength aluminium stabilised conductor is used 
for the superconducting coils to minimise the thickness of 
the support cylinder The cross-sectional view of the con- 
ductor is shown in Fig. [17] A copper matrix NbTi multi- 
filamentary core wire is clad with aluminium stabiliser which 
is mechanically reinforced by means of "micro-alloying" 
and "cold work hardening" [13 14 1 . The overall yield strength 
of the conductor is measured to be higher than 220 MPa at 
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Fig. 16 Gradient magnetic field generated by tlie COBRA magnet. 
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Fig. 17 Layer structure of tlie central coil (left) and the cross-sectional 
view of the conductor (right). Units in mm. 



4.4 K. The measured critical current of the conductor as a 
function of applied magnetic field is shown in Fig. [18] where 
the operating condition of COBRA is also shown (the high- 
est magnetic field of 1.7 T is reached in the coils with the 
operating current of 359.1 A). This indicates that COBRA 
operates at 4.2 K with a safety margin of 40% for the con- 
ductor 



the superconducting magnet is efficiently cancelled around 
the LXe detector since the distribution of the magnetic field 
generated by the compensation coils is similar to that of the 
superconducting magnet around that region but the sign of 
the field is opposite. Figure [19] shows that the expected dis- 
tribution of the residual magnetic field in the vicinity of the 
LXe detector is always below 5 x 10"^ T. 
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Fig. 19 Distribution of the residual magnetic field around the LXe de- 
tector The PMTs of the LXe detector are placed along the trapezoidal 
box shown in this figure. 



3.3 Performance 
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Fig. 18 Measured critical cun'ent of the conductor as a function of 
the applied magnetic field. The operating condition of COBRA is also 
shown. 



3.2.2 Compensation coils 

The COBRA magnet is equipped with a pair of resistive 
compensation coils in order to reduce the stray field from the 
superconducting magnet around the LXe detector (Fig. [15). 
The stray field should be reduced below 5x10"^ T for the op- 
eration of the PMTs of the LXe detector. The stray field from 



3.3.1 Excitation tests 

The superconducting magnet was successfully tested up to 
379 A coil current, which is 5.6% higher than the operating 
current. The temperature and stress distribution are moni- 
tored by using resistive temperature sensors and strain gauges 
attached to the coils, respectively. The temperatures of all 
the coils are found to stay lower than 4.5 K even at the high- 
est coil current. 

The strains of the coils are measured by the strain gauges 
up to 379 A coil current at the central coils and the support 
cylinder where the strongest axial electromagnetic force is 
expected. We observe a linear relation between the strains 
and the square of the coil current, which is proportional to 
the electromagnetic force acting on the coils. This indicates 
that the mechanical strength of the coils and the support 
strength are sufficient up to 379 A. 

The effect of the compensation coils is also measured. 
The gain of the PMTs of the LXe detector, measured using 
LEDs (see Sect. 17.1.11 ). is suppressed by a factor of 50 on 
average in the stray magnetic field from the superconducting 
magnet, whereas it recovers up to 93% of the zero field value 
with the compensation coils. 
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3.3.2 Mapping of magnetic field 

The magnetic field of COBRA was measured prior to the 
experiment start with a commercial three-axis Hall probe 
mounted on a wagon moving three-dimensionally along z, 
r and <p. The ranges of the measurement are |z| < 110 cm 
with 111 steps, 0° < (p < 360° with 12 steps and 0cm < 
r < 29 cm with 17 steps, which mostly cover the positron 
tracking volume. The probe contains three Hall sensors or- 
thogonally aUgned and is aligned on the moving wagon such 
that the three Hall sensors can measure B,, Br and B^ indi- 
vidually. The field measuring machine, the moving wagon 
and the Hall probe are aligned at a precision of a few mrad 
using a laser tracker However, even small misalignments 
could cause a relatively large effect on the secondary compo- 
nents, Br and B^ due to the large main component B-, while 
the measurement of B^ is less sensitive to misalignments. 
Only the measured B. is, therefore, used in the analysis to 
minimise the influence of misalignments and the secondary 
components B, and B^ are computed from the measured B^ 
using the Maxwell equations as 



B^(z,r,^) = B^(za,r,^) + 
Br(z,r,cf>) = Br(zQ,r,(f>) + 



1 r dB,(z' 
r dB,(z' 



d(p 
r,<p) 



-dz' 



dz ■ 



(1) 
(2) 



The computations require the measured values for B,- 
and B^ only at the plane defined by z = Zo- This plane is cho- 
sen at zo = 1 mm near the symmetry plane at the magnet cen- 
tre where Br is measured to be small {\Br\ < 2x10""' T) as ex- 
pected. The effect of the misalignment of the B^-measuring 
sensor on B^{zo, r, (p) is estimated by requiring the recon- 
structed B, and B^ be consistent with all the other Maxwell 
equations. 

Since r = is a singularity of Eqs. ([T]i and (|2|i and the 
effect of the misalignment of the B,-measuring sensor is not 
negligible at large r due to large B,-, the secondary compo- 
nents are computed only for 2 cm < r < 26 cm. The mag- 
netic field calculated with a detailed coil model is used out- 
side this region and also at z < -106 cm and z > 109 cm 
where the magnetic field was not measured. 

The magnetic field is defined at the measuring grid points 
as described above and a continuous magnetic field map to 
be used in the analysis is obtained by interpolating the grid 
points by a B-spline fit. 

4 Drift chamber system 
4.1 Introduction 

The Drift CHamber (DCH) system of the MEG experiment 
ifTSll is designed to ensure precision measurement of posi- 



trons from ju^ — » e^y decays. It must fulfil several stringent 
requirements: cope with a huge number of Michel positrons 
due to a very high muon stopping rate up to 3 x 10^ /i^/s, be 
a low-mass tracker as the momentum resolution is limited 
by multiple Coulomb scattering and in order to minimise 
the accidental y-ray background by positron annihilation- 
in-flight, and finally provide excellent resolution in the mea- 
surement of the radial coordinate as well as in the z coordi- 
nate. 

The DCH system consists of 16 independent modules, 
placed inside the bore of COBRA (see Sect. O aligned in 
a half circle with 10.5° intervals covering the azimuthal re- 
gion € (191.25°, 348.75°) and the radial region between 
19.3 cm and 27.9 cm (see Fig.l20ll. 




Fig. 20 View of the DCH system from the downstream side of the 
MEG detector The muon stopping target is placed in the centre, the 16 
DCH chamber modules are mounted in a half circle. 



4.2 Design of DCH module 

All 16 DCH modules have the same trapezoidal shape (see 
Fig-IZD with base lengths of 40 cm and 104 cm, without any 
supporting structure on the long side (open frame geometry) 
to reduce the amount of material. The modules are mounted 
with the long side in the inner part of the spectrometer (small 
radius) and the short one positioned on the central coil of the 
magnet (large radius). 

Each module consists of two detector planes which are 
operated independently. The two planes are separated by a 
middle cathode made of two foils with a gap of 3 mm. Each 
plane consists of two cathode foils with a gas gap of 7 mm. 
In between the foils there is a wire array containing alter- 
nating anode and potential wires. They are stretched in the 
axial direction and are mounted with a pitch of 4.5 mm. 
The longest wire has a length of 82.8 cm, the shortest of 
37.6 cm. The anode-cathode distance is 3.5 mm. The two 
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Fig. 21 Main components of a DCH module: anode frame with wires 
(front), middle catiiode and hood cathode (back). 



The two detector planes are enclosed by the so-called 
hood cathode. The middle as well as the hood cathodes are 
made of a 12.5 //m-thick poly amide foil with an aluminium 
deposition of 2500 A. 

Thanks to such a low-mass construction and the use of 
a helium-based gas mixture, the average amount of material 
in one DCH module sums up to only 2.6 x 10 "^Xq, which 
totals 2.0 X lO^^Xo along the positron track. 

All frames of the DCH modules are made of carbon 
fibre. Despite being a very light material, it offers the ad- 
vantage of being rather rigid. Consequently, there is only a 
negligible deformation of the frames due to the applied pre- 
tension compensating the mechanical tension of the wires 
(1.7 kg) or the cathode foils (12.0 kg). 



wire arrays in the same module are staggered in radial direc- 
tion by half a drift cell to resolve left-right ambiguities (see 
Fig.l22ll. The isochronous and drift lines in a cell calculated 
with GARFIELD 1 16] are shown in Fig.|23]for B = 1 .60 T. 




Fig. 22 Schematic view of the cell structure of a DCH plane. 



4.3 Charge division and vernier pads 

A preliminary determination of the hit z coordinate, Zwode, 
is based on the principle of charge division. To achieve this 
goal, the anodes are resistive Ni-Cr wires with a resistance 
per unit length of 2.2 ki3/m. Firstly, the z coordinate, is mea- 
sured from the ratio of charges at both ends of the anode wire 
to an accuracy better than 2% of the wire length. 

Secondly the information from the cathodes is used to 
improve the measurement, by using a double wedge or vernier 
pad structure 017lll8lfT9l with pitch A - 5 cm, etched on the 
cathode planes on both sides of the anode wires (see Fig.l24ll 
with a resistance per unit length of 50 i3/m. 



Electron drift lines from a wire 




Fig. 23 Isochrons (dashed green) and drift (solid red) lines of a cell 
withB = 1.60 T. 




Fig. 24 Double wedge or vernier pads etched on the aluminium depo- 
sition of the cathode foil. 



Due to the double wedge structure, the fraction of charge 
induced on each pad depends periodically on the z coordi- 
nate. The period A is solved by means of Zmode- In total there 
are four cathode pad signals for each wire. To increase the 
capability of this method the vernier pads of one cathode 
plane is staggered by AjA in the axial direction with respect 
to its partner plane. 
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4.4 Counting gas 

In order to reduce the amount of material along the positron 
trajectory, the helium-based gas mixture He : CaHe = 50:50 
is adopted as a counting gas ( M20II21II ). For this mixture the 
disadvantages of small primary ionisation and large diffu- 
sion of helium is balanced by the high primary ionisation, 
good quenching properties and high voltage (HV) stability 
of ethane. A reliable HV stability is very important as the 
DCHs are operated in a high-rate environment and with a 
very high gas gain of several times 10^. The main advan- 
tages of this mixture are a rather fast and saturated drift ve- 
locity, ~ 4 cm/fj.s, a Lorentz angle smaller than 8° and a 
long radiation length (Xo/p - 640 m). 

4.5 Pressure regulation system 

A specially designed pressure control system manages the 
gas flows and gas pressures inside the DCH modules and in- 
side the bore of COBRA, which is filled with almost pure 
He to also reduce multiple scattering along the positron path 
outside the DCH volume. The pressure regulation value, de- 
fined as the pressure difference between inside and outside 
the DCH, is on average 1.2 x 10"^ bar regulated to a preci- 
sion better than 0.2 x 10"^ bar. This sensitivity is required 
as larger pressure differences would induce large distortions 
in the thin cathode foils and consequently lead to changes 
in the anode-cathode distances resulting in huge gain inho- 
mogeneities along the anode wire as well as from wire to 
wire. 

4.6 HV system 

The DCHs are operated with a voltage of ~ 1800 V applied 
to the anode wires. In the first years of operation a HV sys- 
tem based on the MSCB system (see Sect. ll0.l~4l ) was used. 
In this case one power supply |[22l was used as the primary 
HV device which fed 16 independent regulator units with 
2 channels each. The primary HV was daisy-chained along 
the 16 regulator modules. 

In 2011 the MSCB. based HV system was replaced by 
a new commercial one. Two plug-in units 1.23 J with 16 in- 
dependent channels each, operated in a Wiener Mpod mini 
crate, were chosen to reduce the high-frequency electronic 
noise contribution (mainly 14 MHz). 

4.7 Readout electronics 

Cathode and anode signals are conducted via printed circuit 
board tracks to the pre-amplifiers, which are mounted at the 
frames of the DCH modules. The cathode signal is directly 



connected to the input of the pre-amplifiers, whereas the an- 
ode signal is decoupled by a 2.7 nF capacitor as the anode 
wire is on positive HV potential. 

The pre-amplifier was custom-designed to match the ge- 
ometrical boundary conditions and electrical requirements 
of the DCH system. The circuit is based on two operational 
amplifiers with feedback loops and protection diodes. The 
total gain of the pre-amplifier is ~50 and the bandwidth is 
140 MHz for the cathode signals and 190 MHz for the anode 
signals. In addition, the anode output is inverted to match the 
required positive input of the DRS read-out chip. The noise 
contribution of the pre-amplifier is 0.74 mV. 

All output signals from the pre-amplifiers are individ- 
ually transferred by non-magnetic, coaxial cables (Radiall 
MIL-C-17/93-RG178) to the back-end electronics. At the 
end-cap of the COBRA spectrometer there is a feedthrough 
patch panel which has resistive dividers splitting the anode 
signals into two outputs in the proportion 1:9. The larger 
signal goes to the DRS input, whereas the smaller one is 
amplified by a feedback amplifier and summed with several 
other anode outputs to form a DCH self-trigger signal. 

4.7.7 Alignment tools and optical Survey 

As the DCH is a position-sensitive detector, the position of 
the wires and the pads must be known very precisely. The 
alignment procedure consists of two parts, a geometrical 
alignment and a software alignment using cosmic rays and 
Michel positrons. 

During the construction of each single frame the posi- 
tions of the anode wires and of the cathode pads are mea- 
sured with respect to an alignment pin located at the bottom 
left edge of each frame. Each cathode hood is equipped with 
two target marks (cross hairs) and two comer cube reflectors 
placed on the most upstream and most downstream upper 
edge of the cathode hood (see Fig. |25]i. After assembly of a 
DCH module the position of these identification marks was 
measured with respect to the alignment pin, which allows 
the alignment of the different frames within the assembled 
DCH module and acts as a reference for the wire as well as 
the pad positions. Even though the chamber geometry is an 
open frame construction a geometrical precision of ~ 50 jum 
was achieved over the full length of the chamber for the po- 
sition of the wires and the pads. 

The support structure between two adjacent DCH mod- 
ules is also equipped with target marks and corner cube re- 
flectors on both the downstream and the upstream sides. 

After the installation of the DCH system inside COBRA, 
an optical survey was performed. The positions of the target 
marks and of the corner cube reflectors on the cathode hoods 
and on the support structure were measured with respect to 
the beam axis and the downstream flange of the magnet us- 
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Fig. 25 Comer cube reflectors and target marks on the downstream 
side of the DCH inodules. 



ing a theodolite (target marks) and a laser tracker system 
(corner cube reflectors) with a precision of 0.2 mm. 

4.7.2 Track-based alignment: the Millipede 

The optical survey measurements described above allow the 
determination of the wire positions in the absolute coordi- 
nate system of the experiment. To achieve a better accuracy, 
we developed a track-based alignment procedure consisting 
of three steps: i) an internal alignment of DCH modules; ii) 
the placement of the so-obtained DCH frame in the spec- 
trometer, with a fine adjustment of the relative DCH-target- 
COBRA position; iii) the alignment of the spectrometer with 
respect to the LXe detector (see Sect. 14.7.4b . 

The first step is accomplished by recording cosmic ray 
tracks during the COBRA shutdown periods, to make the 
internal DCH alignment independent of the magnetic field 
map. Events are triggered by dedicated scintillation Cosmic 
Ray Counters (CRC), located around the magnet cryostat, 
which also provide a reference for the drift-time measure- 
ment. The alignment procedure utilises the reconstructed po- 
sition of hit DCH modules to minimise the residuals with 
respect to straight muon tracks, according to the Millipede 
algorithm |24|. Global parameters, associated with displace- 
ments of each DCH module from its initial position, are ob- 
tained with an accuracy better than 150 /zm for each coordi- 
nate. 

The second step relies on consistency checks based on a 
sample of double-turn Michel positron tracks in the magnet 
volume. The relative DCH-COBRA position is fine-tuned 
until the track parameters related to each single turn (which 
are reconstructed as independent tracks by the fitting algo- 
rithm) are consistent within the fit uncertainties. 

4.7.3 Michel alignment 

The chamber positions from the optical survey are also cross- 
checked by using Michel positrons from pt" decays. Both the 
radial chamber alignment and the alignment in z are anal- 
ysed by looking at the difference between the measured co- 
ordinate and that predicted by a fitted trajectory for each 



chamber plane. Because this diagnostic is insensitive to the 
absolute chamber positions, the overall location of the sup- 
port structure in space is constrained by fixing the average 
radial and z positions over all planes to those of the opti- 
cal survey. Each chamber plane is then shifted to bring the 
mean of its pull distribution closer to the average over all 
planes. The shifts are found to be consistent with the re- 
sults in Sect. 14.7.21 Additional eff'ects are investigated with 
Michel positrons such as tilts in the (r-z) plane, and displace- 
ment in the {x-y) plane. However no such eff'ects have been 
found. 

4. 7.4 Alignment with the LXe detector 

Cosmic rays penetrating both the DCHs and the LXe de- 
tector are used to measure the position of the DCH frame 
relative to the LXe detector. From the distribution of the dif- 
ference between the position of incidence on the LXe detec- 
tor reconstructed by the PMTs and that extrapolated from 
the DCHs, the two detectors can be aligned relative to each 
other with an accuracy of 1 mm in the horizontal and vertical 
directions. 

4.8 Calibrations 
4.&.1 Time calibration 

Systematic time offsets between signals are calibrated by 
looking at the distribution of the drift times. It is a broad dis- 
tribution, whose width is the maximum drift time in the cell, 
but it has a steep edge for particles passing close to the sense 
wire. The position of this edge shows systematic offsets for 
different channels (wires and pads), even in the same cell, 
due to different delays in the electronics. The determination 
of the edge position allows us to realign these offsets. 

An essential ingredient of this procedure is the deter- 
mination of the track time, to be subtracted from the time 
of each hit for determining the drift time. If positron tracks 
from muon decays are used, the time can be provided by the 
track reconstruction itself (see Sect. I4.10t . Yet, a better per- 
formance is obtained if cosmic rays are used, with the track 
time provided with high precision by the CRC and corrected 
for the time of flight from the reference counter to the cell. 
This approach disentangles the calibration from the positron 
tracking procedure, and guarantees a relative time alignment 
at the level of 500 ps, comparable with the counter time res- 
olution. The track reconstruction uncertainty due to these 
residual offset errors is negligible with respect to other con- 
tributions. 

The same procedure can be used for both wire and pad 
signals. Alternatively, to avoid problems due to the smaller 
size of pad signals, wire offsets can be aligned first, and then 
the average time difference of wire and pad signal peaks in 
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the same cell can be used to align the pads with respect to 
the wires. 

4.8.2 z coordinate calibration 

A variety of electronic hardware properties influence the ac- 
curacy of the anode wire and cathode pad measurements of 
the z coordinate. As discussed in Sect. 14.31 the hit z coor- 
dinate is first estimated by charge division using the anode 
wire signals and then refined using the charge distributions 
induced on the cathode pads. The anode charge asymmetry 
is defined as 



node 



mode ' 



(3) 



where Q^J^^'^'' is the measured anode charge at the upstream 
(downstream) end. The cathode pad charge asymmetry is de- 
fined similarly with 2™^)°'''° replacing where cathode 
is either the middle or the hood cathode of the cell. Cal- 
culations show that the cathode charge asymmetry depends 
sinusoidally on z to high precision. One group of calibra- 
tions involves adjusting the parameters entering into the an- 
ode z calculation; this improve the probability of obtaining 
the correct pad period from the initial determination of z. 
The anode z coordinate is given by: 



_ L R 

node ~ \ ^ 1^ 

\2 p 



mode 



(4) 



where L is the wire length, R the pre-amplifier input impedance, 
and p the wire resistivity. 

Figure|26]demonstrates the fundamental calibration tool: 
the measured cathode charge asymmetry versus the mea- 
sured Zmode Coordinate for each cell. The sinusoid period is 
forced to the pad pitch {A - 5 cm) by adjusting the factor 
{k + -) and its phase is forced to be for Zmwde = by cal- 
ibrating the relative gain between anode charges Q^^"''" and 

A calibration of the relative gain between cathode pads 
is performed to achieve optimal z resolution. This proce- 
dure uses plots similar to Fig. |26] except that Zmode is re- 
placed with the more precisely determined coordinate from 
fitted tracks Ztmck- A vertical offset indicates a relative differ- 
ence between upstream and downstream cathode pad gains, 
which is zeroed by adjusting the cathode pad gains. The sys- 
tematic uncertainty in the cathode pad gains from this pro- 
cedure is estimated by comparing the results using Ztmck and 
Zanode as the abscissa, to be crf-^ = 70 fim. 

The chambers, being operated at slight overpressure (see 
Sect. 14. 51 ). experience bowing of the outer cathode hood foils. 
The bulge of the hood is largest near z = and at middle 
cell numbers because the foil is only attached at the edges 
of the chamber This means that the two cathode foils are 
not equidistant from the anode wires, but instead the dis- 
tance varies as a function of z even within the same cell. 
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Fig. 26 A fit to the cathode pad asymmetry versus Zmmde f™ the cath- 
ode hood of a particular drift cell. 



The induced charge, as well as the magnitude of the charge 
asymmetry, also differ between the two planes depending on 
z- 

This effect is calibrated by correcting for the mismatch 
between asymmetry amplitudes. For a given ceU, 
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which is related to the ratio of asymmetry magnitudes for 
geometric reasons, exhibits a clear and expected dependence 
on z, as shown in Fig.|27] The dependence for each wire and 
a preliminary evaluation of the pad z coordinate are used 
to rescale the asymmetries on an event-by-event basis for 
a final evaluation of the pad z coordinate. The charge ratio 
dependence of the charge asymmetry magnitude ratio is ob- 
tained by a scatter plot of those two variables summed over 
all wires and values of z- 




35 40 
|z| [cm] 

Fig. 27 A plot of the charge ratio versus |z| of a particular drift cell. 
This is fit to a quadratic polynomial to correct for the relative cathode- 
to-hood gain. 
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4.9 Tracking 

A chain of software algorithms first measures the arrival 
time and charge on each wire and pad associated with the 
passage of a charged particle through a cell. The arrival times 
and charges are converted into spatial coordinates. The in- 
formation from each cell, referred to as a hit, is cross-checked 
with other hits in the same chamber, and groups of hits con- 
sistent with coming from the same particle are then col- 
lected into clusters. Patterns of clusters consistent with com- 
ing from the same particle are then collected into tracks. 

Finally, a list of clusters associated with each track is 
passed to a Kalman filter to fit a trajectory. The fit incorpo- 
rates average energy loss, multiple scattering, and the mea- 
sured magnetic field. It also attempts to link multiple turns 
consistent with belonging to a single trajectory. The Kalman 
filter provides the track momentum, decay vertex, positron 
emission angles, the predicted impact location at the TC and 
the total path length from the target to the TC. An event-by- 
event indicator of the resolution in each of these quantities 
is also given. 



The result is a radial coordinate resolution of cr^^ — 
210 ^m in the core (87%) and cr^^^ — 780 /jtcl in the tail. 
The design resolution was parametrised by a single Gaussian 
with erf ^ 200 yum. 
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Fig. 28 A diagram of the technique for measuring the intrinsic radial 
coordinate resolution. 



4.10 Performance 

This section reports on the DCH performance regarding spa- 
tial, angular, momentum and vertex resolutions. 

4.10.1 Single Hit Resolution 

We refer to the uncertainty in a single hit position measure- 
ment as the intrinsic resolution, i.e. not affected by multiple 
scattering. A technique for measuring the intrinsic radial co- 
ordinate resolution is illustrated in Fig.|28] The method con- 
siders clusters with exactly one hit in each chamber plane, 
on adjacent cells. A local track circle in the {x,y) plane is cal- 
culated from the position of this cluster and the two neigh- 
bouring clusters. Each hit is then propagated to the central 
chamber plane using this trajectory and the difiference in the 
two radial coordinates at the central plane is interpreted as 
a measure of the radial intrinsic resolution. Figure l29]shows 
the distribution of R plane o - R plane 1 ■ This distribution is fit to 
the convolution of two double Gaussian functions. 

The reason for using a double Gaussian function is that 
the single hit resolution for drift chambers usually has large 
tails, which are due to the non-Gaussian fluctuations of the 
primary ionisation, to the large diff'usion and the irregular 
shape of the time-to-distance relations for hits produced in 
the peripheral regions of drift cells, to wrongly left/right as- 
signments, frequent for hits at small drift distances, to wrong 
drift time measurements in some unusually noisy signals, 
etc. Along with the non-Gaussian material eff'ects (multiple 
scattering and dE/dx), it produces tails in the track parameter 
resolutions (angles, momentum and vertex). 
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Fig. 29 A fit to the distribution of the difference of projected r coordi- 
nates for two-hit clusters. 

A technique for measuring the intrinsic z coordinate res- 
olution is illustrated in Fig. |30] Two-hit clusters belonging 
to a track segment are selected in the same way described 
for the radial coordinate resolution measurement to compare 
the two measured z coordinates. The calculated track circle 
in the {x,y) plane is used to define a local coordinate sys- 
tem, with the centre of the circle defining the origin. Within 
this coordinate system, the polar angle is calculated for 
each overall cluster position and for each hit within the clus- 
ter, as shown in Fig. |30(a)| From the three cluster coordi- 
nates, a quadratic trajectory for z as a function of the local 
(p is computed. This trajectory is used to project the mea- 
sured z coordinates of the two hits in question to a common 
(p and the resulting distribution of the difference of the pro- 
jected z coordinates, ziate - Zeariy, is used to infer the intrin- 
sic z resolution, as pictured in Fig. |30(b)| This distribution, 
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shown in Fig. [ST] is fit with the convolution of two double 
Gaussian functions, resulting in a z coordinate resolution of 
crfCH ^ goo in the core (91%) and crfc« = 2.1 mm 
in the tail. The largest known contribution to the z coor- 
dinate resolution comes from the stochastic fluctuations of 
the baseline in the presence of noise; this is estimated to be 
^'^nohe ~ -'^^ average. The design z coordinate reso- 

lution was erf "-^ - 300 ;um. 



y 




(a) A view of a track segment in (b) A view of a track segment in 
the {x,y) plane. See text for details, the 7/j> plane. See text for details. 



Fig. 30 A diagram of the technique for measuring the intrinsic z coor- 
dinate resolution. 
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Fig. 31 A fit to the distribution of the difference of projected z coordi- 
nates for two-hit clusters. 



According to Monte Carlo studies, this method provides 
a significant overestimate of the true resolution. After cor- 
recting for this, we obtain a single Gaussian resolution crg^ - 
9.4 ± 0.5 mrad and a double Gaussian (pe resolution of cr^^ - 
8.4 ±1.4 mrad in the core (80%) and cr^. = 38 ± 6 mrad 
in the tail, where the errors are dominated by the systematic 
uncertainty of the correction. The Monte Carlo resolutions 
are ag^ ~ 9 mrad cr^ ~ 8 mrad, while the design resolutions 
were crg^^^^ ~ 5 mrad. 

It is also important to stress that these resolutions are af- 
fected by correlations among the other positron observable s, 
which can be treated on an event-by-event basis, so that that 
the effective resolutions determining the experimental sen- 
sitivity are erg. - 8.5 ± 0.5 mrad and cr^^ - 7.7 ±1.4 mrad 
in the core. 

The multiple scattering contribution to these resolution 
is crfl,,0. ~ 6.0 mrad, the rest is due to the single hit resolu- 
tion. 



T3 




50 [mrad] 

Fig. 32 A fit to the distribution of 30, = g' "'"™ - e^'""""' on double- 
turn events. The distribution is fitted with a single Gaussian function 
convolved with itself, and the con'esponding width is shown. 



4.10.2 Angular Resolution 

The resolutions in the measurements of the positron angles 
at the target are measured from data by exploiting events 
where the positron makes two turns in the DCHs. Each turn 
is treated as an independent track, fitted, and propagated to 
the beam line where the track angles are compared. The dis- 
tributions of the difference of the two measured angles in 
double-turn events are shown in Fig.|32]for 6^ and in Fig.l33] 
for The resolution in each turn is assumed to be the same 
and these distributions are fit to the convolution of a single 
(0c) or double (0^) Gaussian functions. These functions rep- 
resent the resolution function of the positron angles. 



4.10.3 Vertex Resolution 

The resolution in the position of the decay vertex on the tar- 
get is dominated by the positron angular resolution. For a 
proper evaluation of the angular resolution function, a pre- 
cise knowledge of the correlations between positron angle 
error and vertex position error is required. The average ver- 
tex position resolutions, however, can be measured directly 
by comparing the projected point of interception at the tar- 
get plane on double-turn events. The difference in vertex z 
coordinates of the two turns is fit to the convolution of a 
single Gaussian function with itself, while that of the ver- 
tex y coordinates is fit to the convolution of a double Gaus- 
sian function with itself. After the Monte Carlo corrections 
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5(|) [mrad] 

Fig. 33 A fit to the distribution of 5<l)e = ^l" - (pl"'' '"™ on double- 
turn events. The distribution is fitted with a double Gaussian function 
convolved with itself, and the coiresponding core and tail widths are 
shown, along with the fraction of events in the core component. 



A complementary approach to determining the positron 
energy resolution is possible by using two-turn events as for 
the angular resolution. Figure |35] shows the distribution of 
the energy difference between the two turns. This is fit to 
the convolution of a double Gaussian function with itself, 
the same shape assumed in the fit of the edge of the Michel 
spectrum. A disadvantage of this technique is its inability 
to detect a global shift in the positron energy scale. This 
technique gives an average resolution of cr^^ = 330 keV 
in the core (79%) and cr^^ = 1.56 MeV in the tail, in reason- 
able agreement with the results obtained from the fit of the 
Michel spectrum. A systematic offset of 108 keV between 
the energies of the two turns also appears; the energy of the 
first turn is systematically larger than the energy of the sec- 
ond turn. A related effect is the dependence of the measured 
Michel edge on 9^. These effects point to errors in the mag- 
netic field mapping. 



are applied, the resolutions are cr,,^ - 1.1+0.1 mm in the 
core (86.7%), cr,,^ = 5.3 + 3.0 mm in the tail and cr,^ - 
2.5 ±1.0 mm. The resolutions for Monte Carlo events are 
o-^'^ = 1.0 ± 0.1 mm in the core and cr^'^ = 2.9 ± 0.3 mm. 
The values of cr,,^ are corrected for the correlation with the 
positron energy assumed to be the signal energy. 

The design resolution was cr,, ~ 1 .0 mm without cor- 
recting for correlation. 

4.10.4 Energy Resolution 

The positron energy resolution is measured with a fit of the 
energy distribution to the unpolarised Michel spectrum mul- 
tiplied by an acceptance function and convolved with a res- 
olution function: 

Probability densityiE'^"""''''') = 
(Michel * Acceptance)(E'J"'') (g) Resolution. (6) 

Functional forms for both the acceptance and the resolu- 
tion functions are based on the guidance provided by Monte 
Carlo simulation. The acceptance function is assumed to be: 



Acceptance(E';"') = ^ ' °" , (7) 

and the resolution function is taken to be a double Gaus- 
sian. The acceptance and the resolution parameters are ex- 
tracted from the fit, as shown in Fig. |34] This gives an av- 
erage resolution of cr£^ = 330 + 16 keV in the core (82%) 
and (Te, - 1.13 + 0.12 MeV in the tail. There is also a 60 
keV systematic underestimation of the energy, to which we 
associate a conservative 25 keV systematic uncertainty from 
Monte Carlo studies. This is to be compared with the reso- 
lution goal of cr^^ = 180 keV (0.8%FWHM). 




Ee (MeV) 

Fig. 34 A fit to the Michel positron energy spectrum. The theoretical 
spectrum (dashed black), the resolution function (dashed blue) and the 
acceptance curve (in the bottom plot) are also shown. 



4.10.5 Chamber detection efficiency 

The relative efficiency of each chamber plane is measured 
as the probability to have a reconstructed hit when its neigh- 
bouring plane in the same chamber has at least one hit as- 
sociated to a track. This probability is called the "hardware" 
efficiency, while the probability to have a hit associated to 
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Fig. 35 A fit to the distribution of dE, = E]" '"™ - El"'' on double- 
turn events. The distribution is fitted with a double Gaussian function 
convolved with itself, and the coiresponding core and tail widths are 
shown, along with the fraction of events in the core component. 



- fast response to be used in the online trigger algorithms 
and to avoid rate effects; 

- online fast and approximate (~ 5 cm) impact point posi- 
tion reconstruction for trigger purposes; 

- capability of determining the positron time with an ac- 
curacy ~ 50 ps 

- good (~ 1 cm) impact point position reconstruction in 
the offline event analysis; 

- reliable operation in a harsh environment: high and non- 
uniform magnetic field, helium atmosphere, possibility 
of aging effects. 

- cover the full acceptance region for signal while match- 
ing the tight mechanical constraints dictated by the DCH 
and COBRA 

These goals were achieved through extensive laboratory 
and beam tests Il25ll26ll27ll28l . 



the track is called the "software" efficiency. Low-efficiency 
planes can be traced to operation at voltages below the nom- 
inal voltage. For the other planes the average hardware ef- 
ficiency is ~ 95% and the average software efficiency is 
~ 90%. 

4.10.6 DCH efficiency 

The absolute efficiency of the DCH system for signal and 
Michel positrons is difficult to measure. Yet it is not used 
in the physics analysis, because the /i+ — > e^y Branching 
Ratio is normalised with respect to the number of recon- 
structed Michel or radiative muon decays, so that the abso- 
lute positron reconstruction efficiency cancels. Nevertheless 
in Fig.[34l the fit to the Michel positron spectrum relative to 
tracks associated to a TC hit returns an estimate of the rela- 
tive efficiency versus the positron momentum. In particular 
the relative efficiency for Michel (48 MeV < Eg) versus sig- 
nal positrons is ~ 0.70. This relative efficiency can also be 
estimated by Monte Carlo simulation to be ~ 0.75 in good 
agreement. 

The absolute efficiency of DCH only for Michel positron 
(48 MeV < £,) as estimated from Monte Carlo is efocH"'""' ~ 
90%, while for signal e^^^^ ~ 83%. The largest source 
of inefficiency for the DCH alone are the tracks emitted at 
6e ~ 90°, where multiple turns generate inefficiency in the 
track-finding algorithm. 

5 Timing Counter 

The Timing Counter (TC) is dedicated to precisely measur- 
ing the impact time of positrons to infer their emission time 
at the decay vertex in the target by correcting for the track 
length obtained from the DCH information. 
The main requirements of the TC are: 



5.1 Concept and design of the Timing Counter 

As mentioned in Sect.[T]and visible in Fig.[Tl the TC matches 
the signal kinematics and is compatible with the mechani- 
cal constraints by the placement of one module (sector) up- 
stream and the other downstream. 

Each sector is barrel-shaped and consists of two sub- 
detectors, mounted on top of each other, with full angular 
coverage for positrons from yU^ — > e^y decays when the y- 
ray points to the LXe detector (see Fig. |36]for a sketch of a 
sector). 

At the outer radius, the longitudinal detector, consisting 
of scintillating bars read out at each end by PMTs, is ded- 
icated mainly to time measurements and, thanks to its seg- 
mentation along (f), also provides a measurement of the posi- 
tron impact (f> coordinate. It also provides a measurement of 
the positron impact z coordinate by exploiting the separate 
time measurements of both PMTs. An approximate estimate 
of these variables is obtained online with fast algorithms for 
trigger purposes. This detector is therefore sufficient to re- 
construct all positron variables required to match a DCH 
track and recover the muon decay time by extrapolating the 
measured time back to the target. 

At the inner radius, the transverse detector, consisting 
of scintillating fibres, is devoted to determining with high 
precision the impact z coordinate to improve the matching 
between the DCH track and TC point. 

Each sector is surrounded by a bag made of a 0.5 mm- 
thick EVAL ® foil that separates the helium-filled volume 
of the magnet bore from the PMTs, which are sensitive to 
helium leakage through the borosilicate window II29II30I . A 
buildup of helium inside the PMTs would lead to internal 
discharge making them unusable. The bag volume is contin- 
uously flushed with nitrogen to guarantee a stable operating 
condition. 
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Fig. 36 Schematic picture of a TC sector. Scintillator bars are read out 
by PMTs. Scintillation fibres are placed on top of the scintillator bars 
and read out by APDs. 



5.2 The longitudinal detector 

Each sector of the longitudinal detector consists of an ar- 
ray of 15 scintillating bars (Bicron BC404 [i31]|25]) with ap- 
proximate square section and size 4.0 x 4.0 x 79.6 cm-'; each 
bar is read out by a couple of fine-mesh 2" PMTs for high 
magnetic fields [32| glued at the ends. The bars are arranged 
in a barrel-like shape to fit the COBRA inner profile with a 
10.5° gap between adjacent bars matching the DCH pitch. 
The scintillator type was chosen from among fast scintilla- 
tor candidates through a series of beam tests by comparing 
their timing performances. The PMTs were chosen to match 
the bar size taking into account the magnetic field strength 
and orientation Il33l as well as the mechanical constraints. 

The number of bars is defined by requirements on the 
size of the trigger regions defined in the (0, 6) plane mapped 
onto the TC; the width corresponds to the gap between 
adjacent bars. 

The bar thickness was chosen to acquire sufficient en- 
ergy deposition from positrons to obtain the design time res- 
olution. 

The bar width is dictated by the requirement of provid- 
ing full efficiency for positrons crossing the TC, leaving no 
gap between bars. Full acceptance of positrons emitted in 
the experiment's angular range defines the bar length. 

Furthermore, the details of the bar geometry are also 
dictated by the mechanical constraints coming from the TC 
transverse detector and the DCH system, resulting in the bar 
geometry in Fig. [37] 

5.2.1 Readout electronics 

The signals from the PMTs are interfaced with the DRS chip 
fast digitiser (see Sect. 19.2.1b with the read-out scheme de- 
picted in Fig. [38] Each PMT signal is passively split into 
three channels with 80%, 10% and 10% amplitude fractions, 
respectively. The highest fraction is sent to a Double Thresh- 
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Fig. 37 The layout of a TC bar-PMT assembly. 



old Discriminator (DTD) specifically designed for time anal- 
ysis purposes. 

DTD is a high-bandwidth, low-noise discriminator with 
two different tunable thresholds: the lower one allows pre- 
cise pulse-time measurements (see Sect. 15.4. 2l l while the 
higher one acts as a veto against low energy hits. The value 
of the high threshold determines the detector efficiency and 
is also related to the trigger threshold as discussed in Sect. 15.41 
When the DTD is fired, a standard NIM waveform is gen- 
erated and then digitised by the DRS boards together with 
a copy of the PMT waveform. The NIM and PMT wave- 
forms are analysed offline to extract the hit time. Several 
algorithms can be applied, the best one (used as default) fits 
the NIM waveform. This method reduces the contribution 
from the DRS jitter and dynamic range to the overall time 
resolution. 

The first 10% copy of each signal is duplicated and fed 
to the trigger system and the DRS: in this stage an active 
splitter described in Sect. I9.2.2l is used for the level trans- 
lation needed to match the input ranges of the trigger (see 
Sect. [HJ and DRS boards. The second 10% signal copy is 
sent to a charge integrator to monitor the PMT ageing. Since 
20 11 we implemented a charge monitor based on pulse anal- 
ysis. Consequently we feed 20% of the PMT signals to the 
DRS/trigger inputs. 
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Fig. 38 Schematic picture of ttie read-out electronics of the TC lon- 
gitudinal detector Signals from each PMT are passively split and fed 
into dedicated electronics channels. Different splitting ratios are used 
to optimise time and pulse reconstruction. The PMT charge monitor 
was deactivated in 201 1 and its fraction is added to DRS/trigger input. 



5.3 The transverse detector 

The transverse detector consists of 256 scintillating multi- 
clad fibres Il34l coupled to Avalanche PhotoDiodes (APD) 
ll35l covering the inner surface of a TC sector 

The small sectional area of the fibres (5x5 mm^ plus 
2x0.5 mm wrapping thickness totalling 6.0 mm in size along 
z) satisfies the required space resolution and matches the 
5x5 mm^ APD sensitive area. In order to comply with the 
mechanical constraints, the fibre ends at the APD side are 
bent in two diff'erent ways: one set has almost straight ter- 
mination while the other is "S" shaped with small curvature 
radii (~ 2 cm) and shows a significant light leakage, recov- 
ered with a high-reflectance wrapping 0361 . 

The APD advantages can be summarised as: small size, 
insensitivity to magnetic field, and fast response. The main 
disadvantage is the low gain, at most ~ 10^ if biased near 
the breakdown voltage, and the high capacitance that make 
them sensitive to noise. For this reason a single APD current 
pulse is read out as voltage across a load resistor and ampli- 
fied. Bundles of 8 APDs are mounted on distinct boards car- 
rying amplifying stages, power supplies, HV bias and ancil- 
lary control signals. Due to mechanical constraints, the fibre 
bundle mounted on a single front-end board is interleaved 
with the fibre bundle connected to the adjacent board: the 
resulting assembly, 16 consecutive fibres read out by two 
adjacent boards, covers 9.6 x 40.0 cm^ and is viewed in the 
trigger system as a single block. 

Each APD has an independent discriminated output, ac- 
quired via a FPGA based board and stored. Therefore the 
transverse detector delivers the list of fired fibres for each 
event. 

The reconstructed coordinate of the hit fibre z/ibre can be 
matched to the Zbar measured by the longitudinal detector 
The distribution of Zbar - z fibre for the downstream sector 
obtained with cosmic rays is shown in Fig. |39] Some tails 



are related to a (9(10%) inefficiency of the APD detector due 
to its geometry: incomplete coverage and 1 mm gap between 
fibres. 
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Fig. 39 Distribution of Zbar ~ z fibre for the downstream sector 



5.4 Commissioning and in-situ performance 

The integration of the TC in MEG requires a set of cali- 
brations to determine the detector parameters optimising its 
performance. Hereinafter only those involving the TC alone 
are described. Additional calibrations, involving other sub- 
detectors (e.g. time alignments with the DCH and LXe de- 
tectors) are described in Sect. l7.3.Tl and in fyf\. 

5.4.1 Gain equalisation 

Due to the non-solenoidal configuration of the magnetic field, 
the inner and outer PMTs are subject to different field in- 
tensities and orientations and therefore operate in diff'erent 
working conditions resulting in different gains. Gains need 
to be equalised for their use in triggering and analysis. 

The PMT gain equalisation procedure exploits cosmic 
rays hitting the TC: due to the uniformity of cosmic ray hit 
distribution along the bars the charge and amplitude spectra 
of the inner and outer PMTs are only marginally affected by 
geometrical effects. It is therefore possible to tune the PMT 
gains acting on the bias high voltages and equalise the peaks 
of the Landau distributions (Fig. l40l i within 15% or better 
In the trigger algorithm a further software equalisation is 
performed online to keep a uniform threshold on the whole 
detector By means of dedicated weights for each PMT a 
percent level uniformity is obtained. 
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Fig. 40 Example of charge spectram acquired with a TC bar, superim- 
posed with a Landau function fit. 




Fig. 41 Graphical representation of the Time Walk effect. 



5.4.2 Time Walk effect 

The Time Walk effect is the dependence on the pulse ampli- 
tude of the threshold-crossing time delay as shown in Fig. 141] 
The analog waveforms sampled at 1 .6 GHz are averaged and 
interpolated separately for each PMT over many events to 
obtain a template waveform, used to evaluate the correction 
required in time reconstruction. In the leading edge region of 
each template waveform, the dependence of the delay time 
versus the pulse amplitude is fitted with the empirical func- 
tion: 



TW(x) ^A + By/x + Clogx 



(8) 



where x is the ratio of the DTD Low Level Threshold (LLT) 
value to the PMT pulse amplitude and A, B, C are the pa- 
rameters fit separately for each PMT. These parameters are 
recalculated every year to account for possible change in the 
PMT pulse shape. 

The correction in Eq. |8]is applied on an event-by-event 
basis. As shown in Fig.|42l this function reproduces well the 
experimental data. 

The optimal LLT value balancing Time Walk immunity 
and noise rejection was found by evaluating the timing res- 
olution of each bar with dedicated tests (see Sect. l5.6l l. 

Figure |43] shows the comparison between the time reso- 
lutions obtained with two different LLT values (10 mV and 
25 mV) on the double bar sample (see Sect. I5.6l l. The dif- 
ference between the time resolutions for the two thresholds 
systematically favours 25 mV. From systematic studies of 
the double bar time resolution for LLT values ranging from 
5 mV to 35 mV in 5 mV step, the optimal value of LLT was 
found to be 25 mV. 

5.4.3 Time offsets 

Inter-bar time offsets due to the electronic chains are evalu- 
ated using cosmic rays. All PMTs are equipped with 10 m 
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Fig. 42 Time delay versus amplitude relation with fit superimposed. 
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Fig. 43 Timing resolutions for different values of low-level threshold 
for double bar events. 



of low-loss signal cable to preserve the leading edge of the 
pulses l38l : the equal length for all PMTs minimises the 
possibility of time offsets drifting during the data-taking pe- 
riod (several months) and, together with a continuous moni- 
toring of the relevant variables, ensures stable operation . For 
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each bar, the distribution of the time difference between the 
inner and outer PMTs is acquired. Due to the uniformity 
and isotropy of cosmic rays this distribution is expected to 
be flat and centred at zero; the mean values of these distribu- 
tions are direct measurements of the relative offsets between 
PMTs on each scintillating bar. The time offsets between 
different bars and between the TC and the LXe detector are 
measured using positrons from Michel decays as well as dif- 
ferent calibration approaches Ii39l . 



of the impact point from Eq [17] with the one obtained from 
EqlBl in Sect. EH It was cr, ~ 7.3 cm in 2009 and has 
improved to cr, ~ 5.0 cm since 2010 after refining the algo- 
rithms. 

From the point of view of efficiency at the analysis level, 
the TC trigger and DTD conditional efficiency studies are 
motivated by the need for all triggered events to fire the DTD 
to allow event reconstruction. The relevant parameter is the 
threshold efficiency eoTD defined, for each bar, as the ratio of 
the number of events with NIM signals to the total number 
of events triggered on the same bar. The high-level threshold 
(HLT) is set slightly lower than the trigger one, in order to 
acquire only events with NIM pulse information. Figure |45] 
presents the efficiency cdtd versus the HLT. Averaging over 
all bars yields £dtd > 99.9% for HLT < 400 mV. 
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Fig. 45 DTD efficiency vs high threshold value. The efficiency plateau 
is reached for HT < 400 mV. 



Fig. 44 Cosmic ray hit maps before and after bar time-oifset subtrac- 
tions. 



Figure |44] shows the effect of accounting for time offsets 
on a two-dimensional (z versus (p) cosmic ray hit map. Using 
the measured time offsets (lower panel), the z reconstruction 
algorithm returns the correct position. 



5.5 Trigger settings 

The TC plays a crucial role in trigger (see Sect. 18.3) : it de- 
tects the positron and provides a preliminary estimation of 
its timing and direction. 

The trigger position resolution is a key parameter for the 
direction match with the y-ray entering the LXe detector and 
is evaluated by comparing the trigger-level reconstruction 



5.6 Timing resolution of TC 

If a positron impinges on a TC bar at time Ttc, the times 
measured by the inner (in) and outer (out) PMTs read as: 



Ttc + bin + TW:,, + ^ 



Veff 
L 



— Ttc + i 



+ TWr.u, + 



- z 



Veff 



(9) 



where bi„ out are time offsets depending on the read-out chain, 
TWin.out are contributions from Time Walk effect, Veff is the 
effective velocity of light in the bar and L is the bar length; 
the z axis points along the main axis of the bar and its origin 
is in the middle of the bar. From Eq.|9]the impact time is : 

TWi„ + TWo,„ L 



Ttc — 



2Veff ' 



(10) 
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A more detailed discussion on the determination of the 
factors b is in Sect l5.4.3l and of the factors TW in Sect l5.4.2l 

The timing resolution is determined from hits on multi- 
ple bars. On a sample of two consecutive bars hit by a posi- 
tron, the time difference AT - T2 - T\ is studied. For each 
bar the impact time is defined as in Eq.[TOl The time resolu- 
tion can be written as 



'^-A/K+^i^^^ (11) 

where ctj 2 are the single bar time resolutions and <j track the 
contribution from the track length spread. Under the hypoth- 
esis that the two bars have similar resolutions cr\ ^ 0-2 and 
neglecting the term due to track propagation, the single bar 
resolution is estimated from cr,/ ^J2. This method overes- 
timates the timing resolution since the track length spread 
term is not corrected for Using triple bar events and evalu- 
ating the quantity: 



AT = T2 



(12) 



the effect of the different path length between bars can be re- 
moved at first order, resulting in a better estimate of timing 
resolution. This is shown in Fig.|46]where the red and black 
markers represent respectively the resolutions obtained in 
triple and double bar samples. However, the triple sample 
has significantly smaller statistics so the double bar sample 
is used as the reference tool for checking the detector per- 
formance. 

A certain degradation of the timing performance between 
the Beam Test (o", ~ 40 - 50 ps) and the Physics Run con- 
figuration (cr, ~ 65 ps) has been observed. A few factors 
contribute to this effect: the need for a lower PMT gain in or- 
der to withstand the high rate and match the dynamic range 
of the DAQ/electronics chain; a slightly higher value for the 
low threshold due to additional noise from the surrounding 
environment; intrinsic uncertainties in both double and triple 
bar estimates; and contributions from DRS calibration. 

5.7 Positron timing resolution 

In the proposal the timing resolution of the positron was as- 
sumed to be dominated by the timing resolution of the TC. 
That was because of the optimistic assumption on tracking 
performances, as well as the neglect of the contributions of 
the calibration of the DRS boards and of TC counter offsets. 
Assuming a small contribution from the uncertainty on the 
track length, the design positron timing resolution was to be 
a,^ = 50 ps. 

In Monte Carlo, which incorporates a fairly precise de- 
scription of DCH and of the positron reconstruction algo- 
rithm, an additional contribution due to track length fluctu- 
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50 ps is present. Added to the TC intrinsic 



ations cr\^ 

resolution crj'^'^'^ - 40 ps, it totals a positron time resolu- 
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Fig. 46 Timing resolution on double (black markers) and triple (red 
markers) bar events, with LLT=25 mV. Due to the different rates for 
two-bar and three-bar events, the sample used here is different than the 
one in Fig. 143 1 



In the data, those contributions are evaluated to be o-f^ - 
60 ps and cr^ = 75 ps. Additionally, an estimated contribu- 
tion from DRS calibration is crf*'^ = 25 ps and from TC 
offset calibration cr™' - 40 ps. The total positron resolution 
is cr,^ = 107 ps. 



5.8 Position resolution 

The positron impact point calculated from the time differ- 
ence between the two pulses is obtained from Eq.|9]as 

z = ^ X (f,„ - f„„, - ibi„ - b„„) - (TWi„ - TWou,)) (13) 

The impact point can also be evaluated using the ratio of 
the charges delivered at the inner and outer PMTs 

e„, ^^G^e^'-e. (14) 

Qou,=EG„,„e'^ (15) 

where E is the energy released inside the bar, Gin^out takes 
into account several contributions (i.e. the scintillator yield, 
PMT quantum efficiency and gain), A^s is the effective at- 
tenuation length of the bar Taking the ratio we obtain: 

Qout Gout 

which leads to: 

Goul 



(16) 



Aeff L Qout 
2 \ " Qin 



In 



(17) 



Moreover, from Eq. [T4]the energy release in the bar can 
be estimated without dependence on z: 



out Gif 



(18) 



tion cr, = 64 ps. 



Note that the combination of Eq.[T7]and Eq. [13] provides 
a way to evaluate the ratio Aefr/vetf for each TC bar Assum- 
ing Veff - 14.0cm/ns, the value of A^s is extracted from 
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Fig. 47 Logarithm of the PMT charge ratio vs the z coordinate mea- 
sured using Eq.^]for bar number 1. The linear fit returns AefF and the 
ratio of PMT gains. 

a linear fit as shown in Fig. |47] The Ass values extracted in 
20 1 1 for all bars fall in the range 40 - 90 cm. This apparently 
large dispersion of values is caused by variation in surface 
reflectivity due to surface machining, hand-made pohshing 
and wrapping, while the bulk attenuation is expected to be 
constant between the bars. 

Two methods for impact point reconstruction are used 
in different stages of the data acquisition chain. The on-line 
algorithm for z reconstruction in the trigger, requiring fast 
response and moderate precision, relies on the charge ra- 
tio method, while the offline analysis, aiming at the best 
possible resolution, exploits the PMT time difference. The 
need for offline calibrations to guarantee the ultimate perfor- 
mance (both for z and time resolutions) is satisfied by using 
several tools, presented in the following. 



5.9 Efficiency 

The TC efficiency for Michel and signal positrons emitted 
from the target can be estimated by Monte Carlo. In the 
acceptance region the efficiency is g^^^"^'^''' ^ 45% for 

Michel (48MeV < positrons and ef^^^^'«™' ~ 59% for 
signal positrons. Such a low efficiency compared to a de- 
sign value > 90%, is due to the DCH support frame that 
intercepts a large fraction of the positrons exiting the DCH 
volume; these interacting positrons may lose a large fraction 
of their energy or change their direction by a large angle due 
to multiple scattering. Therefore they do not reach the TC 



and follow a different trajectory from that predicted by the 
DCH and hence fail the DCH-TC matching cut. 

Using the DCH-only trigger (see Sect. [Q i. the TC ef- 
ficiency versus the reconstructed positron energy was also 
estimated from data, with results compatible with the Monte 
Carlo simulation. 

The total positron efficiency requires a track in the DCH, 
a hit in the TC and that the track matches the TC hit. This 
requirements amount to an efficiency g^*^'*^"^'""' ~ 37% for 
Michel positrons and g^*--'^'*""' ~ 48% for signal positrons, 
as estimated by Monte Carlo. A direct measure from data is 
not possible. 

Absolute detection efficiencies are anyway not used in 
the analysis, as explained in Sect. 14.10.61 therefore the val- 
ues reported are just indicative. 



6 Liquid Xenon detector 

6. 1 Physical requirements vs liquid xenon properties 

Liquid Xenon (LXe), with its high density (and short radia- 
tion length) is an efficient detector medium for y-rays. Both 
ionisation charges and scintillation photons are produced by 
radiation. In MEG, only the scintillation photons are used 
to simplify the detector construction and to utilise a prompt 
response of the detector The general properties of LXe are 
summarised in Table [3] 



Properties 


Value & Unit 


Atomic/Mass number 


54/ 131.293 


Density at 161.4K 


2.978 g/cm' 


Boiling/Melting point 


165.1/161.4 K 


Radiation length 


2.77 cm 


Moliere radius 


4.20 cm 


Scintillation wavelength 


178 nm 


Wpi, for electron / a 


21.6 eV/ 17.9 eV 


Decay time (recombination) 


45 ns 


Decay time (Fast/Slow component) 


4.2 ns/22 ns 


Absorption length 


> 100 cm 



Table 3 Liquid xenon properties: W,,;, is the average energy for scin- 
tillation photon, scintillation wavelength is the wavelength where the 
scintillation spectrum has its maximum. 



The MEG y-ray detector requires excellent position, time 
and energy resolutions to minimise the number of accidental 
coincidences, which are the dominant background process. 
These requirements need high light yield. Fast decay time 
helps to reduce pile-up events. 

On the other hand, the peak of the LXe scintillation emis- 
sion spectrum is in Vacuum Ultra Violet (VUV), A ~178 nm. 
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which requires quartz windows in the PMTs, and the op- 
erational temperature is ~165 K, which requires cryogenic 
equipment. Furthermore, the scintillation light can be eas- 
ily absorbed by impurities like H2O, O2, and N2 etc., which 
must therefore be removed efficiently and continuously. 

6.2 Detector Design 



insulation layer. To reduce the material traversed by inci- 
dent y-rays, the window of the outer vessel consists of a thin 
stainless steel plate (0.4 mm thickness), while that of the in- 
ner vessel is made of aluminium honeycomb panels covered 
with carbon-fibre plates to withstand pressure up to 3 bar. 
The total thickness of the y-ray entrance window is 0.075 
Xq. a picture of the LXe detector is shown in Fig.|49] 



We developed detector prototypes to perform various feasi- 
bility tests. The first used 2.3 ( of LXe with 32 PMTs, and 
the second used 68.6 ^of LXe with 228 PMTs [41 1. The pur- 
pose of the first prototype was a proof-of-principle detection 
of y-rays at ~2 MeV from radioactive sources in LXe with 
adequate resolutions. The second prototype was designed to 
confirm the good resolutions for position, time and energy 
for y-rays at ~50 MeV and to study LXe detector operation 
including the purification system toward a full-size detector. 

Based on the prototype measurements, the final detector 
design employing 900 ( of LXe was completed in 2004. The 
construction began in 2005 and was completed in late 2007. 
A schematic view of the LXe detector is shown in Fig.l48] 



Liquid xenon 




Lateral (Upstream) face 



it 



Lateral (Downstream) face 



Fig. 48 Schematic view of the LXe detector: from the side (left), from 
the top (right). 




Fig. 49 LXe detector under construction. A flange of the outer vacuum 
vessel was not closed yet when this picture was taken. 



The C-shaped structure fits the outer radius of COBRA. 
The active volume is read out by 846 PMTs submerged di- 
rectly in LXe to detect the scintillation light, placed on all six 
faces of the detector (inner, outer, upstream, downstream, 
top, and bottom), with different PMT densities. The active 
volume of the detector is ~ 800 £ and covers 11 % of the 
solid angle viewed from the centre of the stopping target. Its 
depth is 38.5 cm, corresponding to ~ 14 Xq and fully con- 
taining a shower induced by a 52.83 MeV y-ray. 

6.2.7 Cryostat 

The LXe detector consists of the inner and the outer ves- 
sels. The inner vessel holds 900 { of LXe and the PMT sup- 
port structure. The outer vacuum vessel is used as a thermal 



In order to monitor the detector condition, 27 tempera- 
ture, two pressure, and two liquid-level sensors are installed 
inside the detector. To maintain the LXe temperature and to 
enable the recondensation of xenon, a custom-designed 200 
W pulse-tube cryocooler BSI is installed on the top of the 
chimney. It does not produce any electrical noise or mechan- 
ical vibrations, thanks to the cooling principle. In case more 
cooling power is needed, e.g. during liquefaction or purifi- 
cation, liquid nitrogen lines are installed to cool down the 
gaseous volume and the wall of the inner vessel. From the 
other two chimneys, the cables of PMT high voltages, sig- 
nals, and sensors etc. are extracted via feedthroughs. A tur- 
bomolecular pump is connected to each vessel to evacuate 
the detector, and a cryopump and a getter pump are installed 
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in the inner vessel to efficiently remove water vapour, which 
is the most dangerous impurity. 



6.2.2 PMTs 

In total, 846 2" PMTs are internally mounted on all sur- 
faces and submerged in LXe; the support structure is shown 
in Fig.|50] In order to allow the operation in LXe, the PMTs 




Fig. 50 View of the LXe detector inside. A total of 846 PMTs are 
mounted on all surfaces. The face on the left side is the incident one 
for y-rays from the target. 



are designed to withstand up to 3 bar The coverage of the in- 
ner surface with active photo-cathodes is ~ 35%. The PMT 
uses a metal-channel dynode structure to achieve a reason- 
able gain inside a short package (the height is 3.2 cm) and 
a high tolerance for magnetic fields as large as ~ 10"^ T 
B4l . directed either along the transverse or the longitudi- 
nal direction. The PMTs are operated at LXe temperature 
~ 165 K, and are equipped with a quartz window that trans- 
mits VUV photons and a bialkali (K-Cs-Sb) photo-cathode 
sensitive to VUV photons. Aluminium strips are added to 
the surface of the cathode to reduce the sheet resistance, 
which increases at low temperature. Heat dissipation from 
the base circuit is minimised by optimising the resistor chain 
(16 M.Q) with Zener diode protection in the last two dyn- 
odes, which keep the voltage constant under a high counting 
rate. Typical PMTs show Quantum Efficiency (QE) ~ 15% 
and gain ~ 1.8x10'' at 850 V. Inside the cryostat, 3.0 m plus 
1.6 m RG-196A/U coaxial cables are used for the PMTs sig- 
nal, and 3.0 m RG-188A/U coaxial and 1.6 m wires for HV. 
These are connected to the feedthroughs. Outside the cryo- 
stat, signal cables are connected to the splitter boards (see 
Sect. |9.2.2| i, and divided into trigger (see Sect.|8) and DRS4 
boards (see Sect. 19.2. Il l which record all the PMT wave- 
forms. 



6.2.3 Purification System 

As discussed in Sect. 16. II scintillation light may be absorbed 
by impurities in LXe, such as H2O and O2 at ppm level. We 
have developed two purification methods, gas-phase ||4T1 
and liquid-phase [45], to remove those impurities. Gas-phase 
purification removes impurities in Gaseous Xenon (GXe) by 
means of a metal-heated getter A diaphragm pump is used 
to circulate GXe without introducing any additional impu- 
rities. Although the gas-phase purification successfully re- 
duces such impurities, its circulation speed is limited (~ 0.6 
{[h). It turned out during the gas-phase purification study 
that H2O was the dominant component for the absorption. 

Liquid purification was developed to improve the circu- 
lation speed and to remove mainly H2O from LXe by us- 
ing a cryogenic centrifugal fluid pump ||461 and molecular 
sieves (MS 13 A). At the normal working point, the flow rate 
is ~ 70 ^/h. Molecular sieves can absorb more than 24 g of 
water, and the cartridge contains heaters which enable a re- 
generation of the cartridge. 

6.2.4 Storage System 

In addition to the detector, two storage systems for LXe and 
GXe have been developed so that 1000 £ of xenon can be 
stored safely when the detector is not operated [47 1. One 
is a 1000 £ Dewar with a pulse- tube cryocooler as well as 
liquid nitrogen cooling lines and a thermal insulation layer. 
The heat income is estimated to be less than 20 W, and the 
cryogenic tank is designed to tolerate a pressure up to 6 bar. 
As a result, LXe can be stored safely without supplying any 
cooling power for 100 h. The Dewar is connected to the de- 
tector with flexible tubes thermally insulated to allow a rapid 
LXe transfer 

The other storage system consists of eight high-pressure 
gas tanks to store GXe, each of which can contain up to 
360 kg of xenon corresponding to 120 £ of LXe. This can 
be used in long shutdown periods. In Fig. [ST] a schematic 
view of the LXe detector system including purification and 
storage systems is shown. 
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Circulation pump LXc storage tank GXe storage tanks 

Fig. 51 Schematic view of LXe storage and purification systems as 
well as of the detector 
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6.2.5 Detector Operation 

Before starting detector operation, the cryostat is evacuated 
by means of turbomolecular and getter pumps to minimise 
impurities inside the detector Meanwhile, GXe is liquefied 
into the 1000 i storage tank. LXe can be transferred from 
the tank to the detector by means of the pressure diff'erence 
between the pressurised tank and the detector After the liq- 
uid transfer is complete, the pressure of the detector is con- 
trolled to be 1.2-1.4 bar by means of the pulse-tube cry- 
ocooler If the cooling power is too much, heaters surround- 
ing the cryocooler are used to reduce the cooling power, and 
if the cooling power is not sufficient, liquid nitrogen is used 
to sustain the pressure. In 2010, the pulse tube cryocooler 
was sufficient to keep 900 { of LXe in the detector, and ~30 
W heaters are constantly used to produce additional heat. In 
case of a problem, or during a short maintenance period, we 
can transfer the LXe from the detector into the 1000 ( tank 
by means of the centrifugal pump. 

6.3 Reconstruction methods 

The reconstruction methods and algorithms are described. 
In the analysis of the LXe detector, a special coordinate sys- 
tem (u,v,w) is used: u coincides with z in the MEG coordi- 
nate system; v is directed along the negative (/>-direction at 
the radius of the inner face (rin=67.85 cm), which is the di- 
rection along the inner face from bottom to top; w = r - rin, 
measures the depth from the inner face. 

6.3.1 Waveform Analysis 

We extract the pulse time and charge from each waveform. 
A typical waveform is shown in Fig. |52(a)| The baseline is 
estimated by averaging the points in the region before the 
pulse for each channel on an event-by-event basis. The dig- 
ital constant fraction method is used to determine the pulse 
time, defined as the time when the signal reaches a prede- 
fined fraction (30%) of the maximum pulse height. To max- 
imise the signal-to-noise ratio for the determination of the 
charge, a digital high-pass filter, based on a moving-average 
method, is applied as shown in Fig. [52(b)] The charge is es- 
timated by integrating the filtered PMT pulse and is later 
converted into the number of photoelectrons. The number of 
average points is 90-100, corresponding to about a 10 MHz 
cutoff frequency. 

6.3.2 Position 

The y-ray conversion point is reconstructed by using the dis- 
tribution of the light detected by the PMTs near the incident 
position. The three-dimensional position is determined by 




i , , , , i , , , , \ , , , , a 

-600 -500 -40(L. , -300 

Time (nsec) 

(a) 




(b) 

Fig. 52 a) An example of a PMT raw waveform from the waveform 
digitiser. b) An example of a wavefomi with high-pass filter. 

fitting the expected light distribution on the PMTs, calcu- 
lated from the solid angles, to the observed distribution. To 
minimise the effect of shower fluctuations, the fit is done 
twice: the first with ~45 PMTs and the second with ~15 
PMTs around the first fit result. The remaining bias on the 
reconstructed position along u and w coordinates is corrected 
using results from the Monte Carlo simulation. The y-ray 
direction is defined by the line segment between its recon- 
structed conversion point in the LXe detector and the recon- 
structed positron vertex on the target. 

6.3.3 Energy 

The reconstruction of the y-ray energy Ey is based on the 
sum of the number of photons detected by the PMTs cor- 
rected for their geometrical coverage, which depends on the 
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PMT locations: the coverage on the outer face is 2.6 times 
less than that on the inner face. This correction works well 
for almost all events except for very shallow ones (w-y < 
3 cm), which require special treatment because they are very 
sensitive to the relative positions of each PMT and the con- 
version point. In this case, the solid angle subtended by each 
photo-cathode at the conversion point can be a better vari- 
able than the coverage for estimating the collection efficiency. 

In Fig. |53] the linearity of the relation between Ey and 
the number of detected scintillation photons is shown. The 




10 10^ 
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Fig. 53 The number of detected scintillation photons versus the inci- 
dent y-ray energy: 4.4, 11.6 and 17.6 MeV y-rays are measured by 
C-W calibration (see Sect. lT3Tl . while 55 and 83 MeV y-rays by the 
CEX calibration (see Sect. l7.3.2t . Error bars assigned to data points are 
too small to be visible. 

conversion factor from the number of photons to energy in 
order to define an absolute energy scale is obtained from the 
55 MeV y-ray peak taken by the Charge EXchange (CEX) 
caUbration (see Sect. I7.3.2l l. 

6.3.4 Timing 

Each PMT determines the time of arrival of the scintillation 
light from the waveform analysis described in Sect. 16.3.11 
To relate this time to the y-ray conversion time, the prop- 
agation time of the scintillation light must be subtracted as 
well as any hardware-induced time offset (e.g. cable length). 
The scintillation light propagation time is subtracted by ex- 
ploiting the distance between the reconstructed interaction 
point and each PMT assuming an effective light velocity 
~8 cm/ns. After the subtraction of these effects, the y-ray 
conversion time is obtained by combining the timings of 
those PMTs which contain more than 50 photoelectrons and 
calculating the minimum value of defined as a sum of 



the squared difference between each calculated and recon- 
structed time. Typically ~150 PMTs are used to reconstruct 
50 MeV y-rays. The PMTs with large contributions to 
are rejected during this fitting procedure to remove pile-up 
effects. 

Finally, the reconstructed time is corrected by the time- 
of-flight between the positron reconstructed vertex on the 
target and the reconstructed conversion point in the LXe de- 
tector 

6.4 Performance 

The LXe detector performance is described in the following 
and summarised in Table |4] 

6.4.7 Position resolution 

The performance of the position reconstruction is evaluated 
by a Monte Carlo simulation and is validated in dedicated 
CEX runs by placing lead collimators at several positions 
in front of the LXe detector Figure |54] shows the result of 
position reconstruction with a L8 cm thick lead collimator 
installed in front of the LXe detector during the CEX run. 
The lead collimator has two 1 cm wide slits along the u co- 
ordinate at V = cm and v = 6 cm, and this figure shows the 
projection onto the v coordinate. Events with the energy be- 
tween 50 MeV < Ey < %Q MeV, and with the y conversion 
depth Wy > 2 cm are selected. Since the thickness of the lead 
collimator is not sufficient to stop 55 MeV y-rays, the floor 
events which are penetrating the lead collimator are also ob- 
served. A double Gaussian function plus a constant term is 
fitted to extract the position resolution of 54.9MeV y-rays, 
and the results are 6.6 mm, 6.7 mm in this example. The av- 
erage resolution at different positions is 6.9 mm. This result 
contains the effect of the slit width itself and of the spread 
of jfi decay points, and is to be compared with the average 
resolution (6.5 mm) from a Monte Carlo simulation of the 
same configuration. The quadratic difference (L8 mm, ex- 
pected to come from the PMT QEs calibration uncertainty) 
between the data and the Monte Carlo simulation is added 
into the position resolution map built from simulation. 

Taking into account the difference between Monte Carlo 
simulation and the data, the average position resolutions are 
estimated to be cr{u^,Vy) ~5 mm and cr„,^ ~6 mm, respectively, 
comparable with the design position resolutions. They are 
also close to erf"- , ~ 4 mm and cr^'' ~ 6 mm estimated by 
Monte Carlo simulation. 

6.4.2 Energy resolution 

The energy response of y-rays at the signal energy is ex- 
tracted from the CEX calibration. A small correction is ap- 
plied to take into account the different background condi- 
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Fig. 54 Reconstructed position distribution with a lead collimator in 
CEX nans. There are two slits 1 cm wide in the 1.8 cm thick collimator 



tions between the muon and the pion beams, and the opening 
angle between the two y-rays. 

The response function of the detector for monochromatic 
y-rays is asymmetric with a low-energy tail due to mainly 
two reasons. One is the interaction of y-rays in the mate- 
rial in front of the LXe active volume, and the other is the 
shower leakage from the front face. Figure |55] shows the 
LXe detector response to 54.9 MeV y-rays. The distribution 
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Fig. 55 Energy response of the LXe detector to 54.9 MeV y-rays for 
Wy > 2 cm in a restricted range of (uy, Vy). The fitting function is de- 



scribed in the text. The resolution is cte 
4.54%. 



1.56% and FWHMe 



is fitted with an asymmetric function F{x) convolved with 
the pedestal distribution h{x) in the CEX run. F{x) is given 
by 



Fix) 



'Aexp(-^l^) (x>xq + t), 

A exp(-5-(r/2 -{x- xq))) {x < xq + t). 



where A is a scale parameter, xq is a peak position param- 
eter, T is a transition parameter and cte^ is a resolution pa- 
rameter indicating the distribution width on the high-energy 
side. Since F{x) shows the intrinsic resolution of the detector 



without pedestal contribution, it can be used for any realistic 
environment with a different pedestal distribution. 

A 3-dimensional mapping of these parameters is incor- 
porated into the likelihood function for the final analysis 
since they depend on the position of the y-ray conversion, 
mainly on its Wy coordinate. As an example, the average en- 
ergy resolution is measured to be cte^ - 1.6% (3 cm < Wy), 
2.0% (0.8 cm < Wy < 3 cm) and 2.7% (0 cm < < 0.8 cm) 
in 20 11. Except for the acceptance edge along v coordinate 
{cte^ ~ 2.5% with |v| > 68.2 cm and Wy >3 cm), the energy 
resolution depends weakly on the Uy and Vy coordinates. 
This number is to be compared with the energy resolution 
of cr^'- = 1.2% for (2 cm < Wy evaluated by Monte Carlo 
simulation. The reason of this slightly worse resolution is 
not fully understood. The behaviour of PMTs such as gain 
stability, angular dependence etc., or optical properties of 
LXe such as convection might be possible sources. 

The design resolution was cte^ = 1.7% over all Wy. 

6.4.3 Timing resolution 

To investigate the intrinsic time resolution of the LXe detec- 
tor due to photoelectron statistics, two PMT groups (even 
PMT IDs or odd PMT IDs) are defined, and the times of the 
same event are reconstructed by these two groups indepen- 
dently. Then the intrinsic time resolution is estimated by the 
time difference between these two results, being dominated 
by photoelectron statistics, while the effects of electronics, 
the position reconstruction and the event-by-event shower 
spread cancel out resulting in crj^'"^ - 36 ps at 55 MeV. 

The absolute timing resolution is evaluated from the time 
difference between two y-rays, emitted back-to-back from 
the ;r" decay during CEX calibration (see Sect. 17. 3. 2t . recon- 
structed by the LXe detector and by a reference preshower 
counter. Figure |56] shows the measured time difference dis- 
tribution having a resolution of erf = 119 ps, which in- 
eludes not only the LXe detector timing resolution but also 
the contributions due to the uncertainty of the jfi decay po- 
sition and to the reference counter The former is evaluated 
to be cr'l" = 58 ps by the tt" beam spread (~ 8 x 8 mm^ beam 
spot size). The reference counter has two plastic scintillator 
plates and both sides are read out by 2" fine mesh PMTs. The 
timing resolution of each plate is estimated by the time dif- 
ference between the two PMTs, and the timing resolution of 
the counter is estimated by these resolutions by taking into 
account the correlation of the plates. Finally, the resolution 
of the reference counter is evaluated to be cr'^^^ = 81 ps. The 
absolute timing resolution of the LXe detector is estimated 
by subtracting these contributions, resulting in cr,^ = 65 ps at 
55 MeV. Figure|52]shows the energy-dependent timing reso- 
lution of the LXe detector Black squares show the measured 
timing resolution erf , while the red circles show the LXe 
detector timing resolution cr, . The black and the red smooth 
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curves are the fit results, and their functions are shown in the 
figure. A vertical dotted line shows the y-ray signal energy 
(52.83 MeV). The timing resolution improves at higher en- 
ergy, which indicates that the photoelectron statistics still 
contributes significantly. This small energy dependence is 
taken into account in extracting the timing resolution for sig- 
nal y-ray s to obtain crf^'^ - 67 ps. This number is in good 
agreement with the timing resolution of cr^'''^'^ - 69 ps 
evaluated by Monte Carlo simulation. 

The breakdown of the time resolution is as follows: 



phe 



^sho 



(19) 



where ctP'^'' is defined above, crj'''" = 24 ps from electron- 
ics contribution, cr™^ - 20 ps from the y-ray time of flight 
uncertainty (which corresponds to depth reconstruction un- 
certainty), and cr*''" = 46 ps from position reconstruction 
uncertainty and the shower fluctuation. 

The final resolution on the y-ray timing cr,^ is obtained 
combining cr^'' with the additional crj'"' spread due to the 
uncertainty in the muon decay vertex on the target as mea- 
sured by extrapolating the positron at the target plane (see 
Sect. 14.1031 1. This spread is no more than cr',"'' ~ 5 ps giving 
a negligible contribution. The results are cr, - 67 ps and 
o-f^c ^ 69 ps. 

These results are to be compared with the design res- 
olution cr'f^^ ~ 43 ps, that was calculated taking into ac- 
count approximately only the contribution from position re- 
construction uncertainty. Taking into account also the un- 
certainty in the muon decay vertex the design resolution is 
again cr, ~ 43 ps. 
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Fig. 56 Time difference reconstructed by the LXe detector and a ref- 
erence preshower counter for 54.9 MeV y-rays. 
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Fig. 57 Energy-dependent timing resolution of the LXe detector (see 
text for details). 

6.4.4 Detector uniformity and energy scale 

A spatial non-uniformity of the detector, due mainly to its 
geometry, is measured by means of the C-W accelerator 
(see Sect. 17.3.11 ) and CEX (see Sect. 17.3.21) and corrected 
to equalise the detector response. 

The Ey scale is constantly monitored by looking at the 
reconstructed 17.6 MeV energy peak from photons produced 
by C-W protons on Li (see Sect. 17. 3. Ti l. The Ey scale and its 
resolution are cross-checked by fitting the y-ray background 
spectrum to a theoretical spectrum obtained with a Monte 
Carlo simulation including radiative muon decay, positron 
annihilation-in-flight and y-ray pile-up, folded with the de- 
tector resolution (see Sect. I6.4.7t . The energy scale and the 
energy resolution are free parameters in the fit, and we can 
compare the fit results with the calibrated energy scale and 
the measured energy resolution in CEX. The systematic un- 
certainty on the energy scale on the basis of these measure- 
ments is ~0.3%, in which the dominant source is from the 
detector uniformity correction (~0.2%). 

6.4.5 Rejection of pile-up events 

It is important to recognise and unfold pile-up events in a 
high muon intensity environment in order to avoid a loss 
of signal efficiency. At 3 x 10''ju^/s beam rate, around 15% 
of triggered events suffer from pile-up. The y-ray pile-up 
signature is identified by three main methods, the shape of 
the waveform sum of all PMTs, the charge distribution of the 
inner and outer face PMTs, and the x^/NDF distribution in 
the time reconstruction. 
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The first method identifies multiple y-ray events with 
diff'erent timing by searching for peaks in the summed wave- 
form of all PMTs. If a pile-up y-ray is detected by this method, 
its contribution is subtracted by using a template waveform, 
and the waveform of the main y-ray is used for energy esti- 
mation. 

The second method identifies multiple y-ray events in 
different positions by searching the PMT charge distribution 
on the inner and outer faces for spatially separated peaks. 
If the event has two peaks, the main peak is retained after 
removing the secondary y-ray. A pile-up removal method 
was developed, to accomplish this efficiently. A position- 
dependent table containing the average charge of each PMT 
in presence of 17.6 MeV y-rays is prepared beforehand. Once 
a pile-up event is identified, the energy of the event is esti- 
mated by fitting the PMT charges to the table without using 
PMTs around the secondary y-ray. Then, the PMT charges 
around the secondary y-ray are replaced with the charges es- 
timated by the fit. Finally, the energy is reconstructed with 
the replaced PMT charges as described in Sect. 16.3.31 

In the third method multiple y-ray events with different 
timings are identified using the algorithm of the first method. 
If an event has multiple y-rays, the /NDF in the time re- 
construction becomes large. Although the PMTs with large 
contributions to the ^"e discarded when the y-ray tim- 
ing is reconstructed (see Sect. I6.3.4K all the PMTs observ- 
ing more than 50 photoelectrons are used to identify pile- 
up events effectively. If the pile-up is not identified in the 
summed waveform or in the light distribution, and is only 
identified in the timing distribution, the event is discarded. 
Similarly, if the pile-up is only found in the light distribution 
but the subtracted energy is negative or larger than 10% of 
the total energy, the event is discarded as well. This cut has a 
98% efficiency, estimated by the time sideband data around 
the signal region. 

6.4.6 Rejection of cosmic ray events 

Cosmic ray events are rejected using topological cuts be- 
cause these events mostly enter the detector from the outer 
face. Therefore the ratio of the charges collected on the inner 
and the outer faces of the detector is smaller in cosmic ray 
events than in muon-beam events, which enter the detector 
from the inner face. The reconstructed depths of cosmic ray 
events are significantly larger than those of signal y-rays for 
the same reason. Fig. |58] shows the two-dimensional scatter 
plot of the charge ratio vs depth. 

In this two-dimensional parameter space, the selection 
criteria are defined to maximise cosmic ray rejection effi- 
ciency while keeping a signal efficiency of 99%. The cut 
discards 56% of cosmic ray events. 

The combined analysis efficiency of this and the pile-up 
cuts is 97%. 
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Fig. 58 Two-dimensional plot of the charge ratio collected on the in- 
ner/outer faces vs the y-ray interaction depth. Black points show the 
data from dedicated cosmic ray mns, and green ones the signal Monte 
Carlo events. The two lines (blue in the online version) show the selec- 
tion criteria for rejection of cosmic ray events. 



6.4. 7 y-ray background spectra 



o 
a 



B 
s 



10-^ 



10^ 



10 




10' 45 50 55 60 

Energy (MeV) 

Fig. 59 7-ray background spectrum during physics ran in 2010 (see 
text for details). 



In Fig. |59]the y-ray background data are shown, while 
the green line shows the fitted spectrum from the Radiative 
Muon Decay (RMD) plus positron Annihilation-In-Flight 
(AIF) spectrum, the blue line the y-ray RMD + AIF spec- 
trum taking into account the effect of pile-up, the black line 
the cosmic ray spectrum, and the red line the fit result taking 
into account all the background spectra. In the 51-55 MeV 
region with Wy > 2 cm, the RMD + AIF spectrum accounts 
for 94-96% of the total background spectrum, cosmic ray 
spectrum for less than 1 %, and the remaining components 
from pile-up, energy response with higher energy tail, etc. 
for 3-5%. 
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6.4.8 Detection efficiency 

The fiducial volume of the LXe detector is defined as \u\ < 
25 cm, |v| < 71 cm and cm < Wy < 38.5 cm. The y- 
ray detection efficiency is estimated using the Monte Carlo 
simulation, and is confirmed by the measurement in CEX 
runs and by the y-ray background rate in physics runs. The 
low-energy tail in Fig.|55]comes from interactions with the 
material in front of the active volume and shower escape 
from the inner face. On the path of a y-ray entering the 
LXe active volume, there are the drift chamber frame, its 
support structure, the COBRA magnet, an entrance window 
of the LXe cryostat, the PMTs and their support structure, 
etc. The largest source of inefficiency is due to the CO- 
BRA magnet, 14.2%. Within the analysis window 48.0 MeV 
< Ey < 58.0 MeV, the average efficiency in Monte Carlo 
simulation is found to be €y - 65%, taking into account the 
positron event distribution. This efficiency was measured us- 
ing Nal single trigger data in the CEX run. By tagging an 
83 MeV y-ray from decay with the Nal detector, the de- 
tection efficiency of the detector for 55 MeV y-rays is found 
to be = 64 - 67%, consistent with the Monte Carlo esti- 
mation. 

Taking into account the analysis efficiency of 97% as 
discussed in Sect. 16.4.5b . the combined y-ray analysis and 
detection efliciency is ey = 0.65 x 0.97 = (63 ± 3)%. This 
value is to be compared with an expected value ey ~ 60%. 



Energy resolution [m'^ > 3 cm] 1.6% 

Energy resolution [0.8 < Wy < 3 cm] 2.0% 

Energy resolution [0 <Wy < 0.8 cm] 2.7% 

Timing resolution 67 ps 

Position resolution 5(uy, Vy), 6(wy) mm 

Detection efficiency 63% 

Table 4 LXe detector performance summary. 



7 Calibrations 

For each kinematic variable it is necessary to know the reso- 
lution, supplemented by the absolute scale (for the energies) 
or the position of the zero (for the relative time and direc- 
tion) and the degree of stability of these parameters over the 
experiment lifetime. We refer to the knowledge of a single 
quantity at a given time as "calibration", and call the ensem- 
ble of these periodically repeated calibrations "monitoring". 

Continuous monitoring of the apparatus is essential for 
two reasons: firstly for early detection of misbehaviour of 
the apparatus and secondly the knowledge of the parame- 
ters has a direct influence on the systematic uncertainties on 
physical variables, acceptances and thresholds. 
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Fig. 60 Reconstructed position of the 25 Q--sources inside the LXe de- 
tector. In this particular mn the detector was partially filled with liquid 
xenon. The different reconstructed shape of O'-sources in gas and in 
liquid xenon is visible (see text for details). 



The issues of calibration and monitoring were incorpo- 
rated into the design of MEG at an early stage, so that several 
auxiliary devices were designed, assembled and integrated 
for this purpose ||481 . 

The measurement of positron and y-ray kinematic vari- 
ables can be obtained from one detector (the y-ray energy 
and positron momentum) or from the combination of the in- 
formation from several detectors (e.g. the relative time and 
direction). Hence we are prepared to deal with calibration of 
each detector separately, or of multiple detectors simultane- 
ously. 

Firstly the devices involved in the calibration of the LXe 
detector are presented: LEDs and point-like o'-sources inside 
the LXe detector, a movable AmBe source and a 9.0 MeV y- 
ray generator made of a nickel/polyethylene sandwich cou- 
pled to a pulsed neutron generator 

Then the calibration of the spectrometer by means of 
Mott scattering of positrons is described. 

Finally the devices that allow the simultaneous calibra- 
tion of multiple detectors: the 1 MeV Cockcroft- Walton (C- 
W) accelerator and the liquid hydrogen and Nal/BGO setup 
for the Charge EXchange (CEX) reaction are discussed. 

A complete calibration cycle is composed of LED, a- 
source, lithium, boron, neutron generator runs. It takes about 
three hours to complete, and it is usually performed two/three 
times per week. 
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7.1 LXe detector calibration 

Unlike the spectrometer, where the magnitude of the posi- 
tron momentum is determined uniquely by the knowledge 
of the geometry of the DCH and of the magnetic field, the 
value of the y-ray energy must be extracted from the number 
of photons detected by the PMTs surrounding the LXe vol- 
ume, once the proper proportionality factor is known. This 
factor contains the LXe light yield, the photo-cathodic cov- 
erage, the PMT gains and the Quantum Efficiencies (QE). 
All these quantities may depend on time: some more than 
others (e.g. the purity of LXe may change, some PMTs may 
be switched off for some runs). 

For this reason a number of calibration lines are avail- 
able to check the energy scale over the full energy range. 
Low-energy calibrations by means of radioactive sources are 
easier and performed more frequently. Although they are 
of limited use to set the absolute energy scale in the sig- 
nal region, they are helpful in finding gross variations of 
LXe purity. At the opposite end, 54.9 MeV y-rays from tt" 
decays make possible to directly measure the detector re- 
sponse, uniformity and resolution close to the signal energy. 
A drawback of this calibration method is the need to change 
the beam polarity and momentum (from /i^ to n^) and the 
usage of a liquid hydrogen target. 

Table |5]presents a list of these lines, which span a broad 
energy range: 

1 . In the low-energy region 4.4 MeV y-rays from an AmBe 
source and 5.5 MeV a-particles from Am sources de- 
posited on thin wires are used to monitor the PMT QEs 
and the LXe optical properties on a daily basis. In addi- 
tion, 9.0 MeV y-rays from capture by nickel of thermal- 
ized neutrons produced by a neutron generator are also 
available (see Sect. lTTWTSl i. 

2. In the intermediate-energy region a C-W accelerator is 
used, two/three times per week, to accelerate protons, 
in the energy range 400-900 keV, onto a Li2B407 tar- 
get, y-rays of 17.6 MeV energy from ^Li(p,yi-i.b)^Be 
monitor the LXe detector energy scale and resolution, 
while time-coincident 4.4 MeV and 1 1 .6 MeV y-rays 
from "B(/7, y4 4yii 6)'^C are used to intercalibrate the 
relative timing of the LXe detector with the TC (see 
Sect. lTJlT l. 

3. In the high-energy region measurements of y-rays from 
tt" decays produced by the tt" CEX reaction in a liq- 
uid hydrogen target are performed once/twice a year (see 
SectHH. 

7.1.1 LEDs and gain evaluation 

PMT gains are estimated by using blue LEDs immersed in 
the LXe at different positions. In total, 44 LEDs are installed 
in the detector. To minimise the position dependence of the 



detected photons, 1 1 LEDs are flashed simultaneously. In 
dedicated gain measurement runs, performed every second 
day on average, LEDs are flashed at several different inten- 
sities, and the PMT gains are evaluated from the photoelec- 
tron statistics following the method described in f41 1. There 
are two kinds of known PMT gain instabilities, a long-term 
gain decrease, and a rate-dependent gain shift. The former is 
typically 0. 1 %/day during physics runs and 0.4%/day during 
CEX calibration runs (see Sect. l7.3.2l l. 

The latter is observed when starting to use the fi beam, 
the typical shift being ~ +2%. It was found that the LED 
light intensity was stable enough to check the long-term sta- 
bility of each PMT gain. To monitor and correct long-term 
instabilities, constantly flashing LED data have been taken 
during physics runs since 2009 at a rate of ~0.2 Hz: the LED 
peak position for each PMT, and its variation with the time, 
is used to bridge the PMTs gain variation between two con- 
secutive dedicated gain measurement runs. 

7.7.2 Point-like a-sources 

A calibration technique based on a lattice of ^^'Am point- 
like (T-sources was developed and applied, for the first time, 
in a prototype LXe detector ||49l . The sources are prepared 
by fixing small portions of a ^^'Am foil to a gold-plated 
tungsten wire by a thermo-compression method. The result- 
ing wire diameter, after the source mounting, is < 150 fim. 

In the LXe detector we mounted five wires, each one 
hosting five point-like sources of ~ 1 kBq/source, for a total 
activity of a: 25 kBq. The wires are positioned in a staggered 
fashion to optimise the range of angles and distances from 
which they are viewed from the PMTs (see Fig.l60]l. 

In GXe the sources are reconstructed as 3 -dimensional 
spots of ~ 1 cm diameter (the range of a-sources in GXe is 
~ 8 mm). In LXe, due to the much smaller range (~ 40 /vm, 
comparable to the wire diameter) part of the scintillation 
light impinges on the wire itself, which projects a shadow on 
the PMTs opposite to the (T-source emission direction. This 
produces reconstructed positions of (^-sources in the form of 
rings surrounding each wire. It is still possible to identify all 
25 sources and furthermore one can select the preferential 
a-source emission direction (e.g. towards the front face) by 
applying a topological cut on the rings. 

Dedicated a-source calibration runs are coUected on a 
daily basis. Thanks to their different pulse shapes in LXe, 
the trigger is able to classify each event as either low-energy 
y-ray or a-source also during normal data taking (beam-ON). 

The constant monitoring of the combined peak from all 
a-sources is used to control the light yield and transparency 
of LXe. A comparison of the observed light distribution on 
each PMT with the one obtained by a detailed Monte Carlo 
simulation results in a measurement of the LXe attenuation 
length 1,4 IJ . This has turned out to be much larger than the 
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Table 5 Typical calibrations performed to determine the LXe detector performance (energy scale, linearity, etc.) with their energy range and 
frequency. 



Process 


Hiiiergy^ivie v ) 


Frequency 


Charge exchange 


;r p — > ifn 


54.9,82.9 


yearly 




ifi —> yy 






Charge exchange 


n^p —> ny 


129.0 


yearly 


Radiative ju+ decay 


fi^ —> e*yvv 


52.83 endpoint 


weekly 


Proton accelerator 


'LiXp, ri7.6(i4.8))'*Bf; 


14.8, 17.6 


weekly 




"B(p,y4.4rii.6)"C 


4.4, 11.6 


weekly 


Nuclear reaction 




9.0 


daily 


AmBe source 


«Be(Q'3-„A,„«)'2C. 


4.4 


daily 











dimension of our detector during all physics runs (/lAbs > 
300 cm). 

7.1.3 Quantum Efficiency evaluation 

The relative QE of each PMT is also measured by using the 
a-sources. LEDs are not suitable for this task since the LED 
wavelength is different from that of LXe scintillation light, 
and the PMT response depends on the wavelength. 

QI^ is first evaluated from the observed a-peak charge 
compared with that from simulation B9l and then corrected 
by using the y-rays from Li peak events (see Sect. 17.3.1) 
to minimise systematic uncertainties related to the (T-source 
positions and angles between sources and PMTs. The mea- 
sured time dependence of the QE of each PMT has also been 
corrected for. 




E, [MeV] 

Fig. 61 The AmBe y-ray spectmm measured by the LXe detector. A 
Gaussian and exponential fit is superimposed on the green spectrum. A 
Gaussian fit is used after subtracting the exponential natural radioac- 
tivity background. 



7.1.4 Americium/Beryllium source 

Since LXe has a different response for or-particles and y- 
rays, A low-energy calibration point (see Fig. |611 is pro- 
vided by 4.4 MeV y-rays from an Am/Be source. A Gaus- 
sian fit (blue spectrum) is used to extract the peak posi- 
tion and the energy resolution. A 50 kBq AmBe source f50l 
stored in a lead repository can be moved, by means of a com- 
pressed air circuit, to the front of the LXe detector. This al- 
lows a rapid check of the LXe light yield and of the correct 
time constants of the scintillation waveforms of a-particles 
and y-rays in the same energy range. This calibration is pos- 
sible only during beam-OFF periods due to the high rate of 
low-energy y-rays present during normal data taking. In the 
latter condition the method described in the Sect. 17. LSI is 
used. 

^The measured QE is actually the product of the photocathode quan- 
tum efficiency times the collection efficiency of the first dynode. The 
distribution of the measured QE of all PMTs in in the (10-20)% range. 



7.7.5 Neutron Generator 

A 9.0 MeV y-ray line from capture of thermal neutrons in 
nickel ^^Ni(n, yg o)^^Ni by using a Neutron Generator (NG), 
allows a rapid and frequent calibration and monitoring of 
the LXe detector in normal (beam-ON) conditions. The LXe 
detector response can be studied when electromagnetic ra- 
diation associated with the beam (radiative muon decay and 
annihilation in flight) illuminates the detector It is there- 
fore appropriate for checking the stability of the LXe detec- 
tor behaviour in case of variable /v-beam intensity. Neutrons 
of 2.5 MeV kinetic energy are produced by a pulsed D-D 
neutron generator, based on the t/(t/,^He)« nuclear reaction, 
with a tunable frequency (1-100 Hz) and intensity (up to 
2.5 X lO'^n/pulse) lISTI . The neutron generator is placed at 
the centre of a polyethylene block placed in front of the LXe 
detector The calorimeter-facing side consists of alternating 
slabs of nickel and polyethylene. A large fraction of the neu- 
trons are thermalized in polyethylene and then captured by 
nickel nuclei. The sandwich structure is optimised to allow a 
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Nickel 9 MeV y-line 



To LXe 
detector 




Fig. 62 The nickel slabs (red), moderator (blue) and neutron generator 
(green). 



large fraction of the 9.0 MeV y-rays from the capture of ther- 
mal neutrons to reach the LXe detector (see Fig. l62t . This 
structure and the materials are able to accommodate several 
requirements: to minimise the unthermalized neutron flux, 
to avoid background from fast neutron interactions in the 
detector, and to minimise the thermal neutron flux out of the 
moderator-nickel box in order to avoid material activations. 

Nickel was selected instead of other materials because 
of its relatively large probability (~ 34% ll52l ) of emitting a 
single capture y-ray at 9.0 MeV. 

The pulsed mode of the NG is crucial to reach a good 
signal-to-noise ratio when the muon beam is operated. Ther- 
mal neutron capture events follow the NG emission by a typ- 
ical average delay of 50-100 fis. Therefore a special "neu- 
tron generator" trigger is used, in which a 100 ;us gate is 
opened 15 ;us after the NG pulse. In this way the uncor- 
rected y-rays from the beam are reduced, and prompt fast 
neutron interactions are suppressed. The neutron generator 
intensity is tuned in such a way as to have only one photon 
per pulse on average in the LXe detector 

A typical y-ray spectrum is shown in Fig. |63] A double 
Gaussian fit is used to extract the peak position and energy 
resolution. Besides the 9.0 MeV line, the y-ray continuum 
from other Ni level is visible above threshold. Data are ac- 
quired with the muon beam on and off. No energy offset is 
observed between these data sets, if the gain calibration with 
beam on/ofF is applied, proving that the LXe detector is well 
calibrated. Figure |64] shows the consistency of this method 
with the other ones. 



7.2 Monochromatic positron beam 

The abundant positron component of the MEG beam (eight 
times more intense than the yu^ component, but normally 
separated and rejected) is used for selecting a monochro- 
matic positron beam. The positrons hit the MEG target and 
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Fig. 63 The 9.0 MeV 7-ray line from neutron capture in '**Ni as recon- 
structed in LXe detector. 
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Fig. 64 The relative position of the NG 9.0 MeV y-ray line compared 
with other y-ray line relative positions, as a function of time. 



are scattered into the COBRA spectrometer by coherent Mott 
scattering on carbon. 

This method has several potentialities which were exten- 
sively studied with MC simulations: 

- A measurement of the momentum and angular resolu- 
tions of the spectrometer with trajectories (double-turn 
tracks) similar to the ones of the signal. 

- A measurement of the muon polarisation based on the 
comparison between the angular distribution of positrons 
from muon decays at rest and that of the Mott scattered 
positrons. 

- A measurement of the spectrometer acceptance. 

- An independent check of the spectrometer alignment. 

A preliminary experimental study of the positron beam 
properties and of Mott scattering events was performed. The 
average momentum of Mott events is centred close to the in- 
cident beam positron momentum of 53 MeV/c with a mea- 
sured one-sigma spread of 450 keV/c, this being largely due 
to the minimum achievable momentum-slit opening. 

The incident positron beam, associated with prompt pro- 
duction in the proton target, is bunched and synchronous 
with the accelerator 50 MHz RF-signal. This signal can be 
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Fig. 66 The momentum distribution of RF-selected Mott events. 



used to measure the time-of-flight of particles from the pro- 
duction target to the spectrometer and so distinguish beam 
positrons from other backgrounds, as shown in Fig.|65] Here 
a sharp peak due to beam positrons is seen at a relative 
time of -15 ns within the modulo 20 ns timing window 
(50 MHz). The broad peak on the right is associated with 
Michel positrons originating from cloud-muons interactions 
in the MEG target. These muons are produced from pion 
decay-in-flight in the surrounds of the proton production tar- 
get. The smaller peak on the left is consistent with the tim- 
ing expected from pions produced in the production target. 
Finally the constant background is associated with Michel 
positrons originating from muon decays after the Wien fil- 
ter With appropriate timing cuts imposed, the selected Mott 
positron line is shown in Fig. |66] with superimposed a fit 
with two Gaussian, yielding a momentum resolution of 450 
keV/c (one sigma) on the high-energy side. 



7.3 Intercalibration of detectors 

7.3.1 The Cockcroft-Walton accelerator 

A 1 MeV proton C-W accelerator f53\ is coupled to MEG 
for calibrating and monitoring the performance of the LXe 
detector, the relative inter-bar timing of the TC and the rel- 
ative timing between the TC and LXe detector II39II54I . Its 
properties are listed in Table |6] 

The C-W accelerator is placed in a separate and inde- 
pendently radiation-surveyed area, in which it can be ser- 
viced and tested (see Fig.lFTli. 

Calibration measurements are carried out with y-rays 
from two different nuclear reactions: ^Li(p, ■yij,(,)^Be pro- 
duces monochromatic 17.6 MeV y-rays and a broad reso- 
nance centred at 14.8 MeV; "B(/7, 74.4yii.6)'^C produces a 
pair of simultaneous y-rays of 4.4 MeV and 11.6 MeV. The 
two reactions have a threshold at a proton energy of 440 keV 
and 163 keV, with a peak cross section of ~ 6 x 10"^ bam 
and ~ 2 X 10 bam, respectively. To increase the number of 
generated y-rays, especially for the reaction on boron, pro- 
tons of energy higher than the threshold are used. These pro- 
tons are captured during the slow-down process in the target, 
leading to an increase in the event rate without increasing the 
p-heam intensity. 

The proton beam energy is therefore set to 500 keV for 
^L/(/7, ri7.6)^Be and to 900 keV for "B(p,r4.4rii.6)'^C. A 
lithium tethraborate (Li2B407) target (p-target to distinguish 
from the /i-target described in Sect. 12.7) is used to generate 
y-rays from both reactions. 

During calibration runs, the target, contained in a vac- 
uum pipe connecting the C-W accelerator to the MEG area, 
is positioned at the centre of COBRA (see Fig.|62]and |68] |. 
The p-target is oriented at 45° relative to the proton beam 
direction, to reduce the amount of material on the path of 
the y-rays directed to the LXe detector. 
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Table 6 Characteristics of the MEG Cockcroft-Walton accelerator 



Proton beam properties 


MEG C-W 


Energy [keV] 


300-1000 


Energy spread (FWHM)[keV] 


< 0.5 


Angular divergence (FWHM) [mi'ad x mrad] 


< 3x3 


Spot size at 3 meter (FWHM) [cm x cm ] 


< 3 X 3 


Energy setting reproducibility [ % ] 


0.1 


Energy stability (FWHM) [ % ] 


0.1 


Range of the average current [juA] 


1-10 


CuiTent stability [ % ] 


3 


Current reproducibility [ % ] 


10 


Duty cycle [ % ] 


100 




Fig. 68 The Cockcroft-Walton accelerator beam line 



During normal data taking the p-target is positioned down- 
stream outside the COBRA spectrometer. When starting a 
caHbration, the yu-target is removed from the beam Hne by 
means of a compressed hehum system and the p-target is 
inserted to the centre of COBRA by means of an extendable 
bellows system of ~ 2 m stroke. The insertion (or extrac- 
tion) is computer controlled and takes ten minutes. At the 
end of the test the inverse operation is performed, and the ^- 
target reinserted. The reproducibility of its positioning has 
been visually inspected and surveyed to be better than our 
spatial resolutions. 

Steering magnets and monitors are available along the 
proton beam line (see Fig. l68T l for centring the beam on 
the /^-target and for measuring the proton beam properties. 
The data from the lithium reaction are recorded by a low- 
threshold trigger, while a LXe-TC coincidence trigger is used 
to record the two boron y-rays (see Fig,|69]l. 

By means of the C-W calibration lines the energy scale 
of the experiment is constantly monitored, as are possible 
drifts in the relative timing between the LXe detector and 
the TC bars. This allows knowing of the energy scale in the 
LXe detector at a few-per-mil level, and a time alignment 
better than 20 ps. 



7.3.2 n beam and charge exchange set-up 

To calibrate the LXe detector at an energy close to that of 
the signal we use y-rays from neutral pion decay {n^ — > yy). 
A neutral pion is produced in the CEX reaction of negative 
pions on protons at rest n^p — » yf'n. The resulting yf' has 
a momentum of ~ 28 MeV/c in the laboratory frame and 
decays immediately to two y-rays. The photons are emitted 
back-to-back in the tt" rest frame with an energy of 

EX^ — ^ 67.5 MeV. 

^ 2 

In the laboratory frame, the photon energies are 

111 

^r,,2 =r^(l±y6cos0*) (20) 

where /? is the tiP velocity and 9* the center-of-mass angle 
between the photon and the tt" direction. 

Differentiating Eq. (l20l i. the energy spectrum of the two 
photons in the laboratory frame 

dA' dA^ dcosfl* 

= X (21) 

dEy d COS 0* dEy 

is flat (because dN/dcos 0* is constant) between the energies 

EyMin- 2 2 ^|l-p■ ^ ' 

A strong correlation exists between the laboratory relative 
angle between photons and energies. In particular the ex- 
tremal energies are obtained for photons which are emitted 
at relative angle 180° in the laboratory, that is along the 
flying direction. 

This strong correlation between energy and opening an- 
gle can be used to precisely define the energy of one photon 
by tagging the direction of the opposite one. By selecting 
back-to-back y-rays in the laboratory, by means of an auxil- 
iary detector opposite to the LXe detector (an array of Nal or 
BGO in our case), the energies are at the edge of the spec- 
trum, and the higher the collinearity, the more monochro- 
matic the two y-rays. For instance by requiring an opening 
angle larger than 175° (170°) the energy spread cte^IEj is 
better than 0.2% (1%). A few percent of the time one of the 
two y-rays converts into a e^e^ pair, either internally (Dalitz 
decay) or in the target material. 

Additionally the tt" radiative capture on protons pro- 
duces a 129 MeV y-ray coincident with a 9.0 MeV neutron. 

A negative pion beam of 70.5 MeV/c is chosen, partly 
due to the maximum current available for the BTS and partly 
because at this momentum a maximum TOF separation be- 
tween pions and all other particles is achieved at the centre 
of COBRA. The measured intensity is 1 .7 MHz at a proton 
current of 2 mA (lower intensities selectable via slit setting) 
with a round spot-size of 8 mm (sigma), fully illuminating 
the target cell. 
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Fig. 70 Scheme of the hquid-hydrogen target with dedicated transfer 
lines. The 100 litres butfer is filled by means of a standard hydrogen 
bottle up to 2.5 bar: this quantity is sufficient to have a liquid volume 
of 150 cc in the cell. When the cell is cooled down the hydrogen flows 
from the buffer to the target cell. 



A Liquid-Hydrogen (LH2) target is used for the CEX re- 
action. It is a cylindrical cell of 50 mm diameter, 75 mm 
length placed at the end of a 2 m long support. Both the LH2 
target (cold) window and the vacuum (warm) window are 
made of 135 fim Mylar in order to minimise multiple 
scattering and absorption outside the hydrogen-filled vol- 
ume. The LH2 target with dedicated transfer lines is shown 
in Fig. |70] During normal operation the cell is filled with 
150 cc LH2 at 20 K, kept cold by a continuous flow of liquid 
helium. The usage of a mechanical refrigerator is not possi- 
ble due to the high magnetic field. 

The operation of the target is relatively safe, since all 
hydrogen is contained in a ~ 100 i reservoir (the "buff'er") 
at 2.5 bar at room temperature. When liquid helium cools the 
target cell, hydrogen is liquefied at the centre of COBRA. It 
takes three days to set up the LH2 target and a few hours to 
liquefy the hydrogen. A 250 ( liquid-helium Dewar allows 
operation of the target for a little more than 48 hours. 



To select y-rays with a back-to-back topology a sodium 
iodide (Nal) detector (2008-2009) or a Bismuth Germanate 
(BGO) detector (2010-2012) is placed opposite to the LXe 
detector. 

This "opposite side" detector is placed on a movable 
stage (see Fig. 1711 1 that is remotely controlled to allow the 
mapping of the entire front face of the LXe detector When 
displaced from the z = plane, the Nal detector is always 
tilted to face the centre of the LH2 target. 

The Nal detector consists of nine crystals of Nal(Tl) 
scintillator and two layers of lead/plastic scintillator for tim- 
ing measurements (the lead/plastic scintillator can be re- 
moved during energy measurements). Each crystal has a size 
of 62.5 X 62.5 x 305 mm^ (11.8 Xq) and is read out by a 
lOxlOmm^ APD |55J. The lead/scintillator detector is made 
of two plastic scintillators measuring 60 x 60 x 7 mm-' pre- 
ceded by a lead converter 5 mm thick, and is placed in front 
of the central Nal crystal when needed. 

The Nal detector was replaced for the 2010 CEX run 
by a BGO detector composed of 4 x 4 crystals, because of 
the better energy and position resolution of the latter This 
translated into a better monochromaticity in the back-to- 
back topology and into a better efficiency. A typical corre- 
lation spectrum of the energy seen by the LXe detector vs 
the energy seen by the "opposite side" detector in a back-to- 
back topology is shown in Fig.l72l 

7.3.3 Radiative muon decay (RMD) 

Observation of the RMD process — > e^yvv demonstrates 
that the apparatus can detect coincident e-y events, and mea- 
sures directly the global time resolution and the relative off- 
set. At the first stages of the experiment we took runs at low 
beam intensity to reduce the accidental background; later, 
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Fig. 72 Energy recorded by the LXe detector vs the energy recorded 
by the BGO detector during a CEX caHbration run. The extremal 
energies of the calibration y-rays at 54.9 MeV and 82.9 MeV anti- 
correlated in the two detectors are visible. 
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Fig. 73 Distribution of t^y for MEG standard trigger The peak is from 
RMDs, the flat component is from accidental coincidences. 



the time resolution was good enough to also clearly see the 
RMD peak at the nominal beam intensity of ~ 3 x 10^ fi^/s. 
Figure |73] shows the distribution of the e-y relative time fgy 
for RMD decays using the standard MEG trigger 

In order to obtain the resolution on fgy for signals, the 
resolution in Fig|73]must be corrected for the y-ray energy 
dependence shown in Fig|57]and for the positron energy de- 
pendence using Monte Carlo resulting in cr,^^ - 122 ns. 



8.1 Requirements 

An experiment to search for ultra-rare events in a huge beam- 
related background needs a quick and efficient event selec- 
tion, which demands the combined use of high-resolution 
detection techniques and fast front-end, digitising electron- 
ics and trigger The trigger system plays an essential role 
in processing the detector signals to find the signature of 
yt/^ — > e^y events in a high-background environment. 

The most stringent limitation to the trigger latency orig- 
inates from the matching with the DRS (see Sect. 19.2. 11 1, 
whose digitisation needs to be stopped at latest 450 ns af- 
ter the event occurrence. If the trigger were later, the charge 
stored in part of the DRS cells may be overwritten and a few 
samples of the waveform (in particular those recording the 
leading edge of PMT signals) may be lost. The speed ob- 
jective must be achieved while keeping the signal efficiency 
close to unity. On the other hand, the event selection needs to 
provide a significant background rejection. The trigger rate, 
as explained in Sect. 18.51 should be kept below 10 Hz so as 
not to overload the DAQ system, which would result in a 
significant increase in the dead time. 



8.2 The trigger scheme 

In order to accomplish these goals, we developed a system 
consisting of VME boards arranged in a tree-structure (see 
Sect. 18.41 ). based on on-board Field-Programmable Gate Ar- 
rays (FPGA) to process the detector signals and reconstruct 
the observables needed for the event selection li56J . The sig- 
nals are digitised by means of 100 MHz flash- ADCs so as to 
obtain an estimate of the signal amplitude (at the level of a 
few-per-mil accuracy) and timing (with a few ns resolution). 
The advantage of choosing an FPGA-based digital approach 
is manifold: 

- operation of reconstruction algorithms to match with the 
above speed and efficiency requirements; 

- availability of RAM memories for storing the digital data 
stream (for both signal input and algorithm output) into 
cyclic buffers, which makes the trigger system an inde- 
pendent digitiser system as well, to be used as a back-up 
for the DRS; 

- versatility of the trigger scheme, as it is possible to im- 
plement selection criteria, other than for ju^ e^y de- 
cays, related to calibration and other event types; the 
acquisition of pre-scaled triggers associated with those 
events during normal data taking enables us to monitor 
the stability of both detector response and trigger effi- 
ciency as a function of time. 
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8.3 Online algorithms 

The set of observables to be reconstructed at trigger level 
includes: 

- the y-ray energy; 

- the relative e-y direction; 

- the relative e-y timing. 

The stringent limit due to latency prevents use of any in- 
formation from the DCH; the electron drift time toward the 
anode wires may be as long as 200 ns, too long to match 
the latency requirement. Therefore a reconstruction of the 
positron momentum cannot be obtained at the trigger leveS 
Even if there is no explicit requirement on the in the trig- 
ger algorithm, the requirement of a TC hit is equivalent ask- 
ing for a positron with momentum > 45 MeV. 

The y-ray energy is the most important observable to 
be reconstructed. The steep decrease of the spectrum at the 
end-point (equal to the signal energy for ju^ — > e^y) empha- 
sises its role in suppressing the background. For this reason 
the calibration factors for LXe PMT signals (such as PMT 
gains and QEs ) are continuously monitored and periodi- 
cally updated in a dedicated database. The energy deposited 
in LXe is estimated by the linear sum of PMT pulse ampli- 
tudes, which is accomplished by cascaded adder stages (with 
the proper bit resolution) up to the master trigger board. The 
sum of about 250 channels (one for each inner PMT, one 
to four PMTs for other LXe faces) is sensitive to coherent 
noise in the electronics. To cope with that, adder stages are 
preceded by the subtraction of the signal pedestal, updated 
online by the average of out-of-pulse samples. 

The interaction time of the y-ray in LXe is extracted by 
a fit of the leading edge of PMT pulses, with a resolution of 
the order of 2 ns. The same procedure allows us to estimate 
the time of the positron hit on the TC with similar resolu- 
tion. The relative time, a useful tool to suppress our acciden- 
tal background, can thus be obtained from their difference; 
fluctuations due to the time-of-flight of each particle to their 
respective detectors are within the resolutions. 

The amplitudes of the inner-face PMT pulses are also 
sent to comparator stages to extract the index of the PMT 
collecting the most photoelectrons. This provides a robust 
estimator of the interaction vertex of the y-ray in LXe, the 
Moliere radius being similar to the PMT size (see Table |3]l. 
This vertex and the target centre provides an estimate of the 
y-ray direction. On the spectrometer side, lacking any DCH 
information, one has to rely on the coordinates of the TC on- 
line reconstructed hit. Under the assumption of the momen- 

* Although excluded from the e^y tiigger, DCH anode sig- 

nals are collected indeed by the system to provide stand-alone triggers 
associated with tracks in the DCH. For each chamber, anode signals 
are grouped in two sets (as coming out of inner and outer wires) and 
fanned-in linearly before being fed to Typel boards. The track condi- 
tion is met if four out of five consecutive chambers are hit. 



tum being that of the signal and the direction opposite to that 
of the y-ray, by means of Monte Carlo simulated events, the 
PMT index is associated with a region of the TC indexed by 
a bar number and a z-segment. If the online z-coordinate of 
the TC hit (see Sect. I5.5l l falls in this region, the relative e-y 
direction is compatible with the back-to-back condition. 

8.4 Hardware implementation 

The trigger system is organised in a three-layer hierarchic al 
structure, as shown in Fig.|74] which consists of: 

1) VME 6U boards (so-called "Type 1"), each to receive 
and sample 16 input analog detector signals at 100 MHz 
by the 10-bit flash- ADCs |57| to be processed by means 
of digital algorithms implemented on a Virtex XC2 VP20 
FPGAchip|l5l; 

2) VME 9U boards (named "Type 2"), each equipped with 
a Viitex XC2VP20 FPGA chip [58], to collect digital 
pieces of information from Typel boards associated with 
dififerent detectors and combine them to obtain estimates 
of the kinematic observables of interest (namely energy, 
direction and time of flight) of decay particles. 




Fig. 74 Hierarchy of the trigger system: (from left to right) Typel 
boards; Type2 intennediate boards; Type2 master board. 

The input stage of Typel boards consists of plug-in Front- 
End (EE) mini-cards to receive the signals from detectors 
(in AC coupling) via the active splitters (see Sect. |9.2.2] l. A 
low-pass filter with cutoff at 33 MHz (to match Nyquist's 
criterion) is obtained by means of an RC integrator through 
the feed-back branch of an AD8138 differential driver The 
common-mode voltage of the output signal can be further 
adjusted by using a programmable AD5300 DAC to match 
the dynamic range of the FADCs (which is set to 0-2 V in 
the case of LXe inner-face PMTs and 0-1 V for all other 
channels). 
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The FPGA digital output is arranged in a 48 -bit word, 
encoded with charge, time and other pieces of information 
(which depend on the detector) relevant to event reconstruc- 
tion. The output bus is transmitted to the upper-layer boards 
by means of a DS90CR483 Low Voltage Differential Signal 
(LVDS) serialiser clocked at 100 MHz, the same frequency 
as the FADC sampling and FPGA algorithm execution. 

The FPGA inputs and outputs are both recorded into 
cyclic memory buffers, as shown in Fig. |75] The compari- 
son of the output memory content with the result of a sim- 
ulator running on the same input data enabled debugging of 
the execution of reconstruction algorithms. 
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MEMORIES 



ALGORITHM 
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Ancillary interface 



VME Interface 



Fig. 75 FPGA fiiTnware structure. 

At the intermediate layer, Type2 boards collect and com- 
bine the information transmitted by nine Typel boards via 
the companion DS90CR484 LVDS deserialisers. Again al- 
gorithm output arranged in 48-bit buses (up to two per board) 
is transmitted to the trigger master board via LVDS serialis- 
ers and there received by deserialisers. Finally, the task of 
the master is to assemble all event characteristics and gener- 
ate a Stop signal whenever a trigger condition is met. More- 
over, this board is the one to receive Busy signals raised by 
each digitiser during data transfer and to issue a Start com- 
mand as soon as that Busy condition has been released by 
all digitisers. 

The synchronous operation of the system (FADC digiti- 
sation, algorithm execution and data flow) is guaranteed by a 
set of "Ancillary" boards distributing clock signal and con- 
trol (Start and Stop) signals to all digitising electronics. A 
master ancillary board hosts a SARONIX SEL3935 oscilla- 
tor to deliver the 19.44 MHz master clock which, although 
being not essential for trigger synchronisation, is quite im- 
portant for ensuring accurate detector timing. On the trig- 
ger side, each board utilises a PLL (based on a CY7B994V 
Roboclock 1591) to multiply the input clock frequency by a 
factor 5 so as to generate 16 replicas of the 100 MHz clock 
signal throughout a board. The relative phase between these 
can be programmed to compensate for the relative timing 
jitter between chips (and, in the case of the FPGA, even be- 



tween lines in a bus) and ensure proper synchronisation over 
the whole trigger tree. 

Trigger boards are hosted in three VME crates and ac- 
cessed by front-end code running on as many online PCs. 
Each board is equipped with a CoolRunner XC2C384 chip 
[58] as an interface between all the chips and the VME bus. 
Input and output data are read out through a standard, double- 
edge block transfer VME master-slave access to ensure a 
throughput rate close to 80 MB/s. Single word (both read 
and write) VME commands are used to access configuration 
registers to set behavioural parameters (such as thresholds, 
look-up table entries, calibration factors, clock skews, etc.) 
as well as to get the status of control signals (namely the 
Busy) from all the digitisers. 



8.5 Operation and performance 

The y-ray response function is obtained with CEX calibra- 
tion close to the signal energy and continuously monitored 
(although at lower energies) by using C-W proton-induced 
y-ray lines; the resolution is estimated to be 3% at 52.83 MeV 
The equivalent energy threshold for //^ — » e^y events is 
computed by comparing the energy spectrum with the one 
obtained at lower thresholds (from 30 MeV to 40 MeV) 
thanks to a pre-scaled trigger operated in parallel with the 
main one. The ratio of the two spectra is shown in Fig. |76] 
and follows a step behaviour smeared with a Gaussian func- 
tion with an energy resolution in agreement with that esti- 
mated by the fit of the 55 MeV line from CEX data. The 
threshold on the y-ray energy measured online is set to 45 
MeV, which guarantees almost full efficiency (> 99.5%) for 
signal events; on the other side, the single rate drops from 
a few hundred kHz (the bulk of which is related to y-rays 
from radiative muon decays) to 1.5 kHz. 
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Fig. 76 Spectrum of Ey measured offline by the LXe detector, as ob- 
tained with the threshold at 45 MeV on the online measured y-ray en- 
ergy, superimposed to a lower-threshold one. The corresponding effi- 
ciency curve, resulting from the ratio of the two spectra, is also shown. 
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Selection 


Measured (s ' ) 


Expected (s ' ) 




Ey>2 MeV 


3x 10' 


2x 10' 


Ey > 45 MeV 


1.5 X 10^ 


1.3 X 10=* 


At^y < 20 ns 


150 


120 


A9cy ~ n 


12 


10 



Table 7 Trigger rates as a function of progressively applied selection 
criteria (y-ray energy, relative e-y timing and direction) at 3 x 10'/j*/s 
stop rate on target. Measured values agree within 20% with expecta- 
tions based on detector simulation and reconstruction algorithm emu- 
lation. 



To assess the relative e-y timing, a 20 ns-wide window 
is fixed to look for coincidences between the two particles. 
Relative offsets are calibrated by reconstructing events due 
to simultaneous y-rays (one detected in LXe, the other in 
the TC) from ^{p,yy)C (see Sect.|7]i. The resolution is esti- 
mated a posteriori by fitting the edges of the window steps 
with an integral Gaussian function and amounts to 3 ns, in 
agreement with the single-detector timing estimate. The co- 
incidence window is wide enough to achieve full efficiency. 

The e-y Direction Match (DM) enables a further reduc- 
tion of the accidental background trigger rate, which is pro- 
portional to the square of the opening angle [60], by one or- 
der of magnitude at least, as expected by detector and trigger 
simulation, as shown in Table [T] 

In Table [8] a subset of the available selections both for 
physics and calibration are listed following their priority or- 
der. For calibration triggers refer to Sect. 0for details. 

The final goal is to optimise the overall DAQ efficiency, 
defined as the product of the trigger efficiency times the 
DAQ live time. Maximising the former is clearly advisable, 
but may be wasted if accompanied by a significant increase 
in the trigger rate and, as a result, in the dead time. The best 
trade off between the two effects is shown in Fig. |77J due 
to limitations inherent in the DAQ (the digitiser read-out 
lasting ~ 24 ms/event), the TC acceptance region is modi- 
fied until the DAQ efficiency reaches ffie maximum allowed 
value. Therefore the only possibility left to overcome such 
limitations was to reduce the dead time by implementing the 
multi-buffer read-out scheme, as described in Sect. 19.4.11 
After this major improvement, we could loosen the DM- 
match selection criteria and the trigger efficiency increased 
to 97%. 



9 Front-end electronics and Data Acquisition system 

9.1 Requirements 

The Data AcQuisition system (DAQ) is dedicated to record- 
ing the signals from all detectors of the experiment. The re- 



Name 


Ey 


Mteyl 


DM 




MEG 


45 MeV 


10 ns 


narrow 


MEG Low Ey 


40 MeV 


10 ns 


narrow 


MEG Wide DM 


45 MeV 


10 ns 


wide 


MEG Wide Izlte^l 


45 MeV 


20 ns 


narrow 


Radiative Decay 


45 MeV 


10 ns 


no 



Name 



Selection criteria 



LXe alone 
DCH alone 
TC alone 
jfi decay 
a 

LED 

Neutron Generator 
Pedestal 



Threshold on Ey 
DCH hit multiplicity 
TC hit multiplicity 
Coincidence LXe & BOO 
Pulse shape discrimination 
Discrimination on driver 
Discrimination on driver 
Random trigger 



Table 8 List of the most important tiiggers available for physics and 
calibration mns (see text for details). 



Live time - online efficiency plane 
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Online efficiency 

Fig. 77 DAQ efficiency as a function of trigger etiiciency and DAQ 
live time. Superimposed markers correspond to optimal values reached 
during 2008, 2009/2010 and 2011/2012 runs (marked as "Multiple 
Buffer"). 



quire ments for this high-precision and high-rate experiment 
are challenging. Energy has to be measured with an accuracy 
better than a few per mil, even for small signals. The timing 
accuracy of the electronics has to be better than 40 ps, and 
the pile-up detection has to distinguish events separated in 
time down to 10 ns. 

The energy from the LXe detector is determined by adding 
several hundred channels. Even a small level of common 
noise on these channels, such as 50 Hz line noise, would re- 
sult in a huge offset of the summed signal. It has therefore 
been decided to digitise the waveforms of all detectors with 
high resolution and high speed. This allows an efficient pile- 
up detection and an event-by-event baseline measurement to 
eliminate any common low-frequency noise. In order to ob- 
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tain the required timing resolution, the signals must be digi- 
tised with a sampling speed above one Giga-Samples Per 
Second (GSPS). 



lower frequency improves the analog performance of the 
chip. Figure |79] shows the functional block diagram of the 
DRS4 chip. 



9.2 Front-end electronics 

While FADCs with sampling speeds beyond one GSPS are 
commercially available, they are impractical to use because 
of their costs, space and power requirements considering the 
3000 channels of the experiment. An alternative is the use of 
Switched Capacitor Arrays (SCA). They sample input sig- 
nals in analog storage cells at high speed. After a trigger, the 
cells can be read out at much lower speed in the MHz range, 
and digitised by a commercial ADC with typically 12-bit 
resolution. Since no suitable SCA was available when the 
experiment was designed, the design of a new SCA chip ful- 
filling all requirements was started [61 1. The latest version 
of this chip called Domino Ring Sampler 4 (DRS4) has been 
used successfully in the experiment since 2009. 

The sampling frequency was chosen as 1 .6 GHz for the 
TC and LXe detectors, which perform high-precision timing 
measurements, and 0.8 GHz for the DCH, which has less 
stringent timing requirements. 

9.2.1 The DRS4 Switched Capacitor Array chip 

To obtain the high sampling speed of a SCA chip, the write 
switches of the analog sampling cells are operated from the 
inverter chain shown in Fig. [78] and called the domino wave 
circuit. This eliminates the need to use a GHz clock, which 
would be very hard to distribute over thousands of channels. 
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Fig. 78 Simplified schematics of the inverter chain of the DRS4 chip. 

The variable resistors between the inverters, together with 
the parasitic input capacitance of the next inverter, a delay 
element, which is used to vary the sampling speed with a 
control voltage U speed- The analog inputs and outputs of the 
DRS4 chip are differential, which reduces the crosstalk be- 
tween channels. The read-out shift register can be operated 
up to 33 MHz. However, we use only half this frequency, 
since the chip read-out dead time is not an issue and the 




Fig. 79 Functional block diagram of the DRS4 chip. 

The analog control voltage to operate the inverter chain 
is obtained from a Phase-Locked Loop (PLL), which fixes 
the sampling speed with high precision to an external ref- 
erence clock over a wide temperature range. The chip con- 
tains eight data channels and a ninth timing channel used 
to sample the external reference clock directly. Since the 
same domino wave circuit operates all nine channels, any 
timing jitter in that circuit can be measured in the timing 
channel and corrected offline to improve the timing accu- 
racy. The analog output of all nine channels is multiplexed 
into a single external 12-bit ADC, which significantly re- 
duces the cost and space requirements. 

9.2.2 Splitters 

A system of active splitters receives the signals from the 
LXe detector and from the TC through the passive splitters. 
Each splitter board receives 16 signals through coaxial ca- 
bles and makes two copies for the DRS and trigger Differen- 
tial signals are used to improve the noise immunity. Output 
signals are sent to digitization systems through flat cables 
(high-density twisted-pairs cable, 0.635 mm pitch, for the 
DRS and low-density twisted-pairs cable, 1 .270 mm pitch, 
for the trigger). The gain for each channel is one and the 
bandwidth is about 1 GHz for the DRS output and 100 MHz 
for the trigger output. Careful choice of the components and 
their placement ensure a small integral non-linearity. A dedi- 
cated test measurement resulted in a combined non-linearity 
of the splitter system and the DAQ boards of ±0.3%. The 
crosstalk between each channel for the full bandwidth sec- 
tion is <0.8%. Furthermore, the splitter board provides an 
analog sum 4 1 of the LXe detector's lateral face signals, 
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which is used by the trigger system. A block diagram of the 
splitter section is shown in Fig. [80] 
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Test 



DRS Output 
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V 




Sum Output 



Fig. 80 Block diagram of a splitter channel for the LXe detector. 



9.2.3 DAQ boards 

The DAQ boards are based on the VME PMC Carrier (VPC) 
(PCI Mezzanine Card (PMC)) board developed at PSI ||62]| . 
This general-purpose VME board contains two Virtex-II Pro 
FPGAs ||58]| . one of which is used to implement the VME64x 
protocol, while the second contains application-specific code. 
The board contains two PMC slots, which carry two DRS4 
mezzanine boards with a total of 32 channels. The high chan- 
nel density makes it possible to fit 640 channels into a single 
VME crate. The experiment uses three VME crates to host 
1728 channels for the DCH read-out, and two crates to host 
846 channels for the LXe detector read-out and 120 chan- 
nels for the TC read-out. 

9.2.4 Trigger bus 

The synchronisation between the DAQ boards and the trig- 
ger is done via the trigger bus, which is a cable running be- 
tween the master trigger unit and the DAQ boards. The first 
12 lines send signals from the trigger to the boards. They 
contain the main trigger, a synchronisation signal, five bits 
carrying the trigger type and five bits with the least signif- 
icant bits (LSB) of the event number The trigger type is 
one of the 32 types defined in the main trigger (see Sect. |8] 
and Table[8]l. The boards and the front-end program can de- 
cide to handle different events flexibly: for example keeping 
the full waveforms for physics events, while analysing on- 
line the waveforms and storing only integrals and timings 
for calibration events. 

The last eight lines of the trigger bus send signals from 
the DAQ boards to the trigger The main signal is the "busy" 
line, blocking the trigger until the board read-out is finished. 
Since all lines of the trigger bus run between FPGAs, their 
meaning can be changed by modifying the firmware only. 



This has proved to be very flexible and was successfully 
used several times to implement new functionality such as 
a delayed trigger. 

9.2.5 DAQ synchronisation 

While the trigger signal can be successfully used to start the 
read-out of all DAQ boards, the timing it carries is rather 
poor The bus structure of the signal distribution causes this 
signal to arrive at slightly different times at each board. In 
addition the signal is routed through the FPGAs to trigger 
the DRS4 chips. This causes ajitter of the signal of ~ 100 ps, 
which exceeds the experiment timing requirements. To over- 
come this problem, a dedicated high-precision clock distri- 
bution was implemented. A low-jitter 19.44 MHz master 
clock is generated in the trigger and fanned out on a clock 
tree with 200 Low Voltage Differential Signalling (LVDS) 
lines, each of them running to an individual mezzanine board. 
The clock is fed into a clock conditioner chip (63'\ on each 
mezzanine board, which removes most of the jitter picked up 
by the cable, leaving a residual clock jitter below 4 ps. One 
output of the chip is directly digitised by a dedicated timing 
channel in each DRS4 chip and used as a timing reference 
for each event. 

The performance of the system was tested with a 50 
ns square pulse with rise and fall times of 2.5 ns. It was 
passively split and the two copies were sent to two arbi- 
trary channels of the digitising electronics. Their relative 
timing was measured by waveform analysis. In 2009 the 
timing accuracy of the DRS4 chip was found to be ~ 110- 
170 ps. This poor timing had two causes. The first was the 
choice of the external components for the loop filter of the 
PLL (see Sect. 19.2. 11 1. A re-optimisation of these compo- 
nents improved the timing significantly. The second cause 
was the continuous digitisation of the DRS4 chip output by 
the ADC. That generated digital noise which affected the 
timing jitter of the DRS4 chip. A firmware modification was 
then applied to stop the ADC when the DRS4 chip is actively 
sampling data, and starting it only during read-out. Minimis- 
ing these two effects, the timing accuracy became ~ 40-50 
ps, comparable to or smaller than the detector timing resolu- 
tion so as not to degrade the experiment's performance. All 
mezzanine boards were modified in June 2010, so all data 
since then were taken with good timing resolution. 



9.3 DAQ System 

The data acquisition system (see FiglSTli is dedicated to read- 
ing out the front-end electronics, combining the information 
into events and storing them on disk. It consists of a cluster 
of PCs connected to the front-end electronics and an event- 
building PC connected through a Gigabit Ethernet switch. 



47 



Signal inputs 

i 



Signal inputs 

Ul iJUiU 




Fig. 81 Simplified schematic 
the lower one just outside the 



of the DAQ system. The individual components are described in the text. The upper part resides next to the detector, 
expeiimental area. The electronic logbook, the system messages and the alann system are omitted for clarity. 



They run the DAQ software, which handles all aspects of 
event-based data acquisition and slow control. 

9.3.1 Online cluster 

The VME boards of the front-end electronics and trigger 
system are located in nine VME crates. Since they produce 
a significant amount of data, each crate is read out by a sep- 
arate PC via the Struck SIS3 100/1 100 VME controller over 
a 15 m optical cable |64|. The PCs contain two dual-core 
Xeon CPUs running at 2.8 GHz under Scientific Linux 5.1. 
The multi-core functionality is important since the PCs that 
read out the DRS chips have a high computing load for the 
online calibration and zero suppression of the raw waveform 
data. If the front-end program run on a single core, the max- 
imum event rate would be only 12 Hz. Running separate 
threads for the VME readout, calibration and network trans- 
fer of calibrated events increases the rate to 38 Hz. The PCs 
have redundant dual power supplies and hot-pluggable mir- 
rored hard disks, so that the failure of a power supply or 



a disk would still allow data taking to continue. The front- 
end PCs are connected to the central back-end PC via a Gi- 
gabit Ethernet switch. This PC uses the same technology, 
but it has more disks. A Redundant Array of Independent 
Disks (RAID 5) system with a net capacity of 2 TB can 
buffer about 2-3 days of data taking. A dedicated process 
copies the data from these disks constantly to the central PS I 
computing centre, where a dedicated computer cluster anal- 
yses the data (offline cluster). The offline cluster contains 16 
nodes, each with four CPU cores, and a disk system with a 
total of 280 TB, and is used both for data analysis and for 
simulation. 

9.3.2 MIDAS software 

The data acquisition software used for the MEG experiment 
is based on the Maximum Integrated Data Acquisition Sys- 
tem (MIDAS) |65|. MIDAS consists of a library written in 
C and C-n-, several tools and a general program framework 
to assist the development of experiment-specific front-end 
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programs as well as analysis programs for online and offline 
use. The MIDAS library implements the following function- 
ality: 

- System Layer. All Operating System (OS) specific func- 
tionalities are covered by a set of functions using con- 
ditional compilation. In this way the upper level of all 
MIDAS programs can use these general functions and 
compile on different OSes. Over the last twenty years, 
many OSes have been supported. Some of them have 
been phased out such as MS-DOS and VMS, others such 
as OSX (Darwin) have been added. Since the system 
layer is relatively small (about 7% of the MIDAS code), 
new systems can be added relatively quickly. 

- Online Database (ODB). The online database is a hi- 
erarchical database which stores all relevant information 
in a central location. For performance reasons it is im- 
plemented completely in shared memory, so that all ex- 
periment processes can exchange data very rapidly at 
an access rate ~ 1 MHz. The MEG ODB is about 5 
MB in size and shared by about 20 processes. The ODB 
contains configuration parameters for the trigger system, 
DAQ VME boards and event building, as well as the cur- 
rent values for the slow control data and event statistics. 
A special feature of the ODB are "hot-links". A process 
can register to be notified through a call-back routine 
when a specific value in the ODB changes. 

- Buffer Manager (BM). The buffer manager is a set of 
library functions, over which processes can exchange 
event-based data. Buffers are segments of shared mem- 
ory to which processes can attach as data producers or 
consumers. Data is passed very efficiently in a first-in 
first-out (FIFO) manner either locally or over the net- 
work using a dedicated TCP/IP protocol. On the Gigabit 
Ethernet network of the MEG experiment a net transfer 
rate of 98 MB/s has been measured. Each consumer can 
either register to receive only a subset of all events (such 
as for monitoring) or request all events (such as the log- 
ger writing data to disk). In the latter case the consumer 
automatically slows down the data producer in case the 
data cannot be digested as fast as produced (commonly 
known as "back-pressure"). 

- System Messages (SM). Using the BM with a dedicated 
buffer, each process sends or receives text messages to 
inform the users about errors or status changes such as 
run start/stops. All SM are centrally recorded during the 
duration of the experiment. 

- Slow Control System (SC). MIDAS contains a frame- 
work to access "slow control" devices (see Sect. [TOl i. 
The framework takes care of transferring information 
between a device and the ODB (read-out) and between 
the ODB and the device (control) by means of a device 
driver implementing two functions for reading and writ- 
ing. 



- History System (HS). While the ODB contains only the 
current value of all slow control values, it is desirable 
to see changes of the variables over time. The history 
system implements a database for storing and retrieving 
all slow control variables over the duration of the experi- 
ment. An optimised binary database is capable of writing 
up to 10^ variables per second to disk, and retrieving the 
time variation of any variable over days or months in a 
couple of seconds. 

- Event Builder (EB). The event builder is a dedicated 
process which combines event fragments from nine dif- 
ferent front-end computers into a single event. Software 
and hardware counters are used to ensure that only frag- 
ments belonging together are combined. In case of a 
fragment mismatch the event builder stops the current 
run and issues an error message. 

- Logger. The logger is an experiment-independent pro- 
gram which stores the events produced by the EB on 
disk. In addition, it writes the slow control data from the 
ODB into the HS utilising the above described hot-link 
scheme. The logger contains basic run control function- 
ality such as stopping a run after a predefined number of 
events and automatically starting the next run. The log- 
ger writes data to a local disk array. A second logging 
process copies the data through a FTP interface to a cen- 
tral computing cluster, from where they are transferred 
to a tape archive. 

- Alarm System. The alarm system is a part of the MI- 
DAS library which runs in each process attached to the 
experiment. It checks individual values in the ODB or 
for the presence of a process. If a value exceeds a pre- 
defined limit or if a process dies, the alarm system can 
produce a system message, stop the run and execute a 
shell script, which in turn can send emails or text mes- 
sages to mobile phones. 

- Online Analyzer. An analysis program based on the 
ROME framework II66II67I runs in the online environ- 
ment sampling and analysing a subset of the data stream. 
It combines a full physics reconstruction with a single 
event display. Using various panels, sub-detector hit pat- 
terns as well as short-term histograms can be inspected 
by the shift crew to identify potential problems. A sec- 
ond copy of the program runs on the offline cluster to 
analyse all events in near-to-real time. Histograms con- 
taining all events are available 10-20 minutes after the 
data has been taken, and allow monitoring of effects which 
become visible only with adequate statistics. More infor- 
mation about the analysis software can be found in 1681 

ii- 

- Web Interface. A dedicated web server (mhttpd) is used 
to control the whole experiment through web browsers. 
It contains an interface to the ODB, so the experimenter 
can view or change any ODB variable. In addition, the 
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interface implements panels for run control (start/stop a 
run plus run status display), system messages, the alarm 
system, the history system and a process interface to start 
and stop individual programs. While the mhttpd program 
is experiment independent, "custom pages" can be de- 
fined to monitor and control a specific part of the ex- 
periment (see Sect. 110.2.2b . The web interface is also 
available from off-site locations for remote monitoring. 

- Run control. While the standard MIDAS system has ba- 
sic run control functionality, more sophisticated run se- 
quencing is needed. For example, for the LED calibra- 
tion of the LXe detector one needs a sequence of several 
runs, each with a diff'erent LED flashing amplitude. To 
fulfil these needs, a dedicated utility called "runsubmit" 
was written. This tool reads an Extensible Markup Lan- 
guage (XML) script and executes a sequence of tasks. It 
can start and stop runs after a certain number of events, 
and change values in the ODB. This allows the tool to 
change the trigger settings and prescaling factors, con- 
figure the data read-out mode and even change slow con- 
trol values such as high voltages or beam line settings. 
Using the sequencing, several runs can be taken auto- 
matically with different settings. 

- Electronic Logbook (ELOG). Instead of using paper 
logbooks, the collaboration uses an electronic logbook 
based on the ELOG program |70|. This program, orig- 
inally an integral part of MIDAS, is now available as 
a stand-alone package. The usage of an electronic log- 
book has proved extremely helpful for communication 
between the on-site shift crew and ofi'-site collaboration 
members, as well as for the retrieval of old logbook en- 
tries tagged with certain keywords. An additional feature 
is a shift checklist, filled out once by each shift, requir- 
ing the inspection of critical elements in the experimen- 
tal area not covered by automatic sensors. This checklist 
is a web form integrated into the ELOG program and 
filled out using a tablet computer in the experimental 
area. This allows a truly paperless operation of the ex- 
periment. 

The MIDAS software runs very stably; in MEG the DAQ 
system is only stopped for configuration changes. If left in 
production mode, it runs for weeks, making most efi'ective 
use of the beam time. 



9.4 Performance 

The performance of the DAQ system can be defined as the 
fraction of recorded events for a given trigger rate. This in- 
cludes dead time caused by the read-out, lost events dur- 
ing data transport, events not recorded between runs, and 
events lost during the downtime of the DAQ system due to 



re-configurations or crashes. Muon beam downtime and cal- 
ibration runs do not contribute to DAQ inefliciency. 

9.4.1 Dead time during event read-out 

The DRS is stopped on each trigger and read out before 
being re-activated. One waveform contains 1024 samples, 
each stored in a 16-bit value, digitised at 16.5 MHz. Each 
DRS chip uses four channels for detector waveforms and 
one for the global reference clock. Two DRS chips are read 
out through a multiplexer by a single ADC, which requires 
61 ns X 1024 x 5 x 2 = 625 ^ls. Each VME board contains 
two mezzanine boards each having four DRS chips. This 
amounts to 1024 (bins) x 5 (channels/chip) x 4 (chips/mez- 
zanine) X 2 (Bytes/sample) = 40 kB per event and mezza- 
nine board. A full VME crate contains 20 boards with 40 
mezzanine boards, which totals to 1.56 MB per crate and 
event. The SIS3 100/1 100 VME interface has a setup time 
of 125 yus per DMA transfer and then a line speed of 83 
MB/s. The transfer of the 40 data blocks therefore takes 
125yusx40-i- 1.56MB/83MB/S = 23.8 ms. The data read- 
out rate from the trigger system is smaller and therefore doe s 
not pose a bottleneck. Before 20 1 1 , the front-end electronics 
was only re-enabled after all data were read out through the 
VME bus, thus causing a dead time of 23.8 ms + 0.625 ms - 
24.4 ms per event. This limited the maximum event rate 
to 41 Hz and caused a live time of 83.2% at our average 
event rate of 6 Hz, given the Poissonian time structure of 
event triggers. At the end of 2010, a multi-buffer read-out 
scheme was deployed in the firmware of the DRS and trigger 
DAQ boards. This scheme uses a FIFO buffer between the 
front-end read-out and the VME transfer, capable of hold- 
ing up to three events for the DRS part and up to four events 
for the trigger part. This allows us to re-enable the front- 
end electronics after the read-out time of 625 and do the 
data transfer in the background. If another trigger happens 
during the VME transfer of the data, it is written into the 
FIFO buff'er without interfering with the data currently be- 
ing transferred. The probability that three or more events 
occur during the 23.8 ms of the VME transfer is 5 x 10 '* at 
a trigger rate of 6 Hz and can be neglected. The dead time 
is therefore determined by the probability to have another 
event during the DRS read-out time of 623 ;us, which causes 
a dead-time of 0.4% and leads to a live time of 99.6%. 

9.4.2 Other sources of dead time 

The data read out through VME are transferred through Gi- 
gabit Ethernet to the back-end PC. Since the TCP/IP proto- 
col is used, no packets are lost. Since a network switch with 
enough buffers and backplane bandwidth has been chosen, 
the probability of data collisions and necessary network re- 
transfers is very small. Constant data rates of 98 MB/s over 
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hours at the input port of the back-end computer have been 
observed, which are more than enough for the typical over- 
all data rate of 20 MB/s. At the back-end PC, the event frag- 
ments are combined into full events. The average event size 
of a typical run is 2.4 MB, consisting of ADC waveforms 
from the trigger system (12%), DRS waveforms for the LXe 
detector (30%), the TC (8%) and the DCH (45%). The re- 
maining 5% covers special waveform channels for timing 
synchronization and the slow control data such as high volt- 
age currents, temperatures and so on. The full event is then 
compressed using the bzip2 algorithm, which results in a fi- 
nal event size of 0.9 MB (62% compression ratio). A special 
parallel version of this algorithm is utilized that makes use 
of all eight CPU cores, so that the full 10 Hz event rate can 
be compressed in real time. 

A typical MEG run takes about six minutes. This time 
was chosen to produce raw data files ~ 4.8 GB large. The 
time it takes to stop and start a new run is six seconds, lead- 
ing to an additional dead time of 1.7%. Although this time is 
small, there are plans to decrease it further by using multi- 
threaded control programs for starting and stopping runs. 



10 Slow Control 

Slow Control in this context means the monitoring and con- 
trol of the physical quantities in the experiment, that vary 
"slowly" over time, such as temperatures, pressures or high 
voltages. It is typically sufficient to read these values once 
every few seconds and record them together with the data 
from the various detectors. This allows us to find during the 
offiine analysis any correlation between changes in detec- 
tor performance and variations in ambient quantities such as 
the temperature. Other quantities need adjustments during 
the data taking: the beam line needs to be tuned differently 
for calibrations and for data taking and the high voltages for 
DCH and PMTs need to be changed from time to time in 
order to compensate for drifts in gain and efficiency. 

Most experiments use many different systems for the 
slow control subsystems. This approach has the problem of 
combining the data from different subsystems into the cen- 
tral DAQ system. While this is viable in large collaborations, 
it is difficult in smaller experiments such as MEG to main- 
tain many different subsystems. Thus the collaboration de- 
cided since the beginning to design a common slow con- 
trol system for the whole MEG experiment. The challenge 
of such a design is to cover all requirements of the exper- 
iment with a single approach. A field bus was clearly nec- 
essary to connect spatially separated sensors, and local in- 
telligence was required for control loops such as automatic 
high-voltage ramping and trip recovery. The system had to 
be flexible and easily extendable, very reliable and still op- 
erational with one part under maintenance. 



10.1 MSCB System 

The new system was named Midas Slow Control Bus (MSCB), 
indicating the close integration into the MIDAS DAQ sys- 
tem. It consists of a communication layer through an RS- 
485 field bus and through Ethernet, and uses 8-bit micro- 
controllers as communication end points and for local con- 
trol loops. A crate-based high-voltage system was designed 
based on this concept. A general-purpose control box SCS- 
2000 was built which can be equipped with various daughter 
cards for up to 64 channels of various analog or digital in- 
puts/outputs. 



10.1.1 Communication 

Physical quantities in MEG are measured at different lo- 
cations. To collect them in a centralised system, an elec- 
tronic connection between the sensors and the DAQ is neces- 
sary. Since analog connections are prone to electronic noise 
pickup, it is desirable to digitise the quantities close to the 
sensors and transfer the information digitally. For the MSCB 
system we choose a field bus based on the RS-485 stan- 
dard. This bus allows a multi-drop topology, where many 
devices are connected in parallel to the same physical bus, 
which uses a 5 V TTL differential signal. The used baud rate 
is 115.2 kBaud/s, which is low enough so that bus lengths 
below 100 meters do not require termination. Measurement 
nodes were designed to use a 8-bit microcontroller which 
implements the MSCB protocol and connects to the serial 
bus via a RS-485 transceiver The low-level bus protocol 
follows exactly the RS-232 standard (one start bit, one stop 
bit, no parity), so that the internal microcontroller UART 
(Universal Asynchronous Receiver/Transmitter) is used. On 
the other side the microcontroller communicates to various 
ADCs and DACs through different protocols such as SPI 
(Serial Peripheral Interface bus) or I2C (Inter-Integrated Cir- 
cuit), as shown in Fig.[82l 
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Fig. 82 Simplified communication schematics of the MSCB system 
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The used RS-485 bus transceivers (MAX1483) have a 
reduced bus load, so that up to 256 nodes can be connected 
to a single physical bus. To connect to the data acquisition 
PC, a so-called "submaster" was built to act as a master on 
that bus segment. On one side it connects to the MSCB bus, 
on the other side it encapsulates all MSCB telegrams into 
network packets and communicates with any PC through 
standard lOBase-T Ethernet. The Ethernet submaster con- 
cept has proved to be extremely versatile. MSCB nodes be- 
longing together logically are connected to one MSCB bus 
segment and interfaced via a dedicated Ethernet submas- 
ter The MEG experiment uses a total of 15 submasters, 
which are connected to the various subsystems and to the 
central experiment network switch. If one segment is down 
for maintenance, the other segments can still be operational, 
which was very useful during commissioning phases. The 
submasters can also be accessed remotely for debugging, 
development and safety, which can happen in parallel to the 
normal DAQ operation. In case the central DAQ system is 
down, a remote backup PC can temporarily take over the 
slow control monitoring, without the need to change any 
physical cable or reboot any system. 

10.1.2 Communication protocol 

The protocol for the MSCB system was designed for min- 
imal overhead. A sequence of just a few bytes is enough 
to read or write a slow control value. Therefore it is possi- 
ble to monitor many hundred nodes per second in spite of 
the relatively slow baud rate of 1 15.2 kBaud on the RS-485 
bus. Each packet starts with a one-byte command, followed 
by optional arguments and a 8-bit CRC code. Usual com- 
mands are reading and writing of values, obtaining informa- 
tion about variables, rebooting a node or updating the node 
firmware through the RS-485 bus. 

Since many nodes can be connected to the same bus, an 
addressing scheme has to be used to communicate with in- 
dividual nodes. Each node contains a unique 16-bit address 
stored in the EEPROM of the node microcontroller To ad- 
dress a certain node, an address packet containing the node 
address has to be sent over the bus. If a node recognises its 
own address, it switches to the "addressed mode" and inter- 
prets all consecutive commands until another node becomes 
addressed. 

Figure |83] shows an example of the communication be- 
tween a submaster and a node. A first packet addresses the 
node, while a second packet requests a certain slow control 
value from that node. The reply to this command contains a 
command code indicating an acknowledge, followed by the 
actual value and the CRC code. 

The MSCB protocol is encapsulated by the Ethernet sub- 
master into UDP packets. A software library on the PC side 
communicates with the submaster using these UDP packets. 
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Fig. 83 Typical communication protocol example between a submas- 
ter and a node. 



The library is written in C and available for common op- 
erating systems such as Linux, Windows and OSX. In addi- 
tion, virtual instruments (VI) have been written to export the 
functionality to the Lab VIEW system from National Instru- 
ments. All MSCB commands require an acknowledge. The 
commands writing a value get a one-byte acknowledge from 
the node after applying the new value, while the read com- 
mand returns the actual value read. This scheme protects not 
only lost or corrupted packets on the RS-485 bus, but also 
lost UDP packets. If an acknowledge is not returned in a 
certain time window, the packet is re-sent several times with 
an increasing timeout. This protocol has been optimised for 
local area network connections, where timeouts of the order 
of a few hundred milliseconds are needed to achieve maxi- 
mal throughput. Tests have shown that this is more effective 
than a TCP connection. A simple password authentication 
scheme ensures that none of the several experiments at PSI 
using the MSCB system talk to each other accidentally. The 
connectionless protocol ensures that the system will recover 
automatically after a network failure or a reboot of a node or 
the submaster 

The MSCB system acquires data and controls values by 
exchanging certain values called network variables. Each 
MSCB node can implement up to 256 network variables. 
Each variable has a name, a type (such as byte or 4-byte 
floating point) and a physical unit. This variable definition 
resides on each node and can be queried using the MSCB 
protocol. Since all necessary information can be obtained 
directly from the MSCB nodes, it is possible to use general- 
purpose programs to connect to the MSCB system and dis- 
play all network variables of a certain submaster without us- 
ing a central configuration database. Adding new nodes or 
removing nodes is trivial. The MIDAS web interface (mhttpd) 
described in Sect. 19.3.2! contains a dedicated MSCB page 
which can be used as an user interface to all nodes of the 
MEG experiment, as in Fig. [84] 
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Fig. 84 Web interface page to the MSCB system. The left column 
shows all available submasters, the second column all nodes connected 
to the selected submaster, and the right column all network variables 
of the selected node. 



10.1.3 Local intelligence 

Each MSCB node and submaster requires a microcontroller 
to implement the MSCB protocol and to realise local intel- 
ligence needed e.g. in control loops. The requirements for 
those microcontrollers are low power consumption, small 
size, integrated analog electronics such as ADCs and DACs 
and user reprogrammability. This lead to the C808 1F12 1 mi- 
crocontroller ITTiI used in most MSCB nodes, which con- 
tains 128kB of EEPROM, 8kB of RAM, a 12-bit ADC and 
two 12-bit DACs and runs at a clock speed of 100 MHz. 

The microcontroller is programmed in C using the Keil 
PK51 development kit |72]. The software is divided into a 
framework implementing the MSCB protocol and general 
housekeeping functions and the application-specific part. Both 
parts communicate via network variables described above. 
The available program space of 128kB is more than enough 
to accommodate the C runtime library including floating 
point functions, the MSCB framework consisting of about 
3500 lines of C code, and still leaves space for several thou- 
sand lines for application-specific functionality such as reg- 
ulation loops. The DCH pressure control system for exam- 
ple uses elaborate code including sensor linearisation, alarm 
functions and Proportional-Integral-Differential (FID) loops 
to keep the gas pressure stable to better than 10'^ bar, all 
implemented in about two thousand lines of C code. 

One special feature of the MSCB system is that it can be 
reprogrammed over the network. The microcontrollers can 
write to their own code space from the user application. A 
small boot loader has been realised and is located in the up- 
per memory of the code space, implementing a subset of 
the MSCB protocol. Using this boot loader, the main pro- 
gram space can be re-programmed. Firmware updates of all 
MSCB devices are therefore very easy. Since the protocol 
is encapsulated fully in the Ethernet communication layer. 



the update can even be sent through the Ethernet submaster 
from remote locations. A typical update contains 50 kB of 
code and requires ten seconds. After the update, the node 
reboots in about 0.5 seconds, reducing downtime to a mini- 
mum. 



10.1.4 High-voltage system 

The MEG experiment requires high-voltage supplies for the 
PMTs of the LXe and the TC, the APDs of the TC and the 
DCH. The MEG experiment requires an accuracy and long- 
term stability better than 0. 1 V, which was hard to achieve 
with commercial systems at the time of the design of the 
experiment. Therefore a high-voltage supply based on the 
MSCB system was designed. Since the experiment has many 
channels which need about the same high voltage, it was 
easier to make a voltage regulator which is fed from a com- 
mon external source. All 848 channels for the LXe detector 
are for example powered from a single high-current com- 
mercial high-voltage source, and the voltage is adjusted for 
each channel separately by a MSCB voltage regulator The 
regulator uses a string of optocouplers, each capable of sus- 
taining 240 V. A total of twelve optocouplers allows voltage 
regulation from zero to 2800 V. The LEDs of the optocou- 
plers are controlled by a regulation loop which compares the 
output voltage with the demand value obtained from a DAC, 
as shown in Fig. [85] 

The output voltage is measured by means of an ADC and 
a voltage divider consisting of two resistors with 41 MOhm 
and 33 kOhm. The current is determined by measuring the 
voltage drop across a resistor Rsense with a diff'erential ampli- 
fier. Using 24-bit sigma-delta ADCs (AD7718 from Analog 
Devices), very good resolutions can be achieved. The ADCs 
and the DAC are controlled by the microcontroller through 
the SPI bus via optocouplers for noise immunity. The ac- 
curacy depends only on the resistors of the voltage divider 
and the current sense resistor, which can be calibrated using 
an external multimeter The calibration is then stored in the 
nonvolatile memory of the microcontroller Accuracies of 30 
mV in voltage and a few nA in current have been achieved 
at 2000 V output voltage. 

10.1.5 RS-232 and GPIB interface nodes 

Although the MSCB system handles many sensors and ac- 
tuators directly, some devices, e.g. COBRA, have their own 
control systems. These devices typically have a RS-232 or 
GPIB interface. To integrate these devices into the MEG 
slow control system, MSCB interface nodes were designed 
for the RS-232 and the GPIP busses, containing only the 
microcontroller with the RS485 interface chip (ADM2486 
from Analog Devices), and a second interface chip for the 
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Fig. 85 Simplified schematics of the high-voltage regulator. The 240 
V Zener diodes protect the optocouplers, so that no more than 240 V 
drop across the transistors of the optocouplers. 

RS-232 bus (MAX232 from Maxim) or the GPIB bus (di- 
rectly connected to the CPU). A 5V DC-DC converter gal- 
vanically decouples the bus power from the local power 
This allows the nodes to be powered from the bus without 
the danger of introducing ground loops between different de- 
vices through their earthing system. 

The interface nodes have two modes of operation. The 
so-called transparent mode encapsulates the device commu- 
nication into MSCB data packages. The PC side contains 
routines which allow direct writing or reading from the RS- 
232 or GPIB bus. The advantage of this mode is that new 
devices can be connected and tested without modification 
of the MSCB firmware since this system works only as the 
transport layer The disadvantage is that each device needs a 
dedicated driver on the PC side which understand the device- 
specific ASCII protocol. 

The second mode encapsulates the device specific pro- 
tocol in the firmware of the interface node. If connected to a 
power supply for example, the voltage and current of the 
power supply are polled periodically and put into MSCB 
network variables. These variables can then be read out with- 
out the knowledge of the device-specific protocol to query 
the voltage and current of the power supply. Parameters in 
attached devices can be changed in a similar way. Writing 
to the associated network variable triggers a callback rou- 
tine in the microcontroller firmware which generates a write 



command to the device through the RS-232 or GPIB bus. 
Using this mode, all devices attached to the MSCB system 
can be talked to in the same way (namely network variables) 
independent of the specific device. This mode has therefore 
been chosen for all devices in the MEG experiment. 

70.7.6 SCS-2000 control box 

The MSCB-based SCS-2000 control box was designed for 
analog and digital inputs and outputs. The box contains an 
Eurocard-size main board with a microcontroller, a battery - 
buffered real-time clock and a 1 MB static RAM. User in- 
teraction happens through a LCD display and four push but- 
tons. The SCS-2000 unit can either be inserted in a 3 HE 
crate or be housed in a stand-alone box with an external 24V 
power supply. A single box can accommodate eight daugh- 
ter cards, each with eight input/output channels of various 
kinds. Analog input cards typically host a 24-bit sigma-delta 
ADC with eight input channels. Other cards were developed 
for analog and digital input and output, including a capaci- 
tance meter for the liquid xenon level gauge or a lOV LED 
pulser for the LXe detector calibration. 

A generic firmware scans all plugged-in daughter cards 
and displays their parameters on the LCD display. Using the 
four buttons on the front panel, the user can scroll through 
all values and change values of those cards which support 
digital or analog output. This is very useful for controlling 
devices locally with an attached SCS-2000 without the need 
of using a laptop or other computers. The parameters are 
visible as network variables and are accessible through the 
MSCB bus. 

10.2 DAQ Integration 

Most slow control parameters in the MEG experiment are 
acquired locally by SCS-2000 units or other MSCB nodes, 
but they need to be integrated into the full DAQ system. 
Measured values have to be written to disk together with the 
event data. Furthermore, the temporal changes of parame- 
ters have to be accessible by the shift operators in order to 
ensure the stability of the experiment. This is achieved by 
a dedicated slow control front-end program, together with 
the MIDAS online database and history system, which are 
described in Sect. 19.3.21 

10.2.1 MIDAS slow control front-end 

A dedicated process runs on the back-end computer which 
connects to all slow control submasters through the Ether- 
net network. All controlled devices are partitioned into sep- 
arate entities called equipments, such as high-voltage sup- 
plies, beam line magnets or ambient temperature sensors. 
Each equipment is read out periodically with a dedicated 
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thread inside the front-end process for optimal read-out fre- 
quency. The readout threads copy all parameters into the 
ODB, where they become immediately available to other 
processes or to the shift operator via the mhttpd web inter- 
face. Control parameters are placed into the ODB as well, 
and linked to their devices via hot-links (see Sect. |9.3.2| l. 
If the user or a process changes a control parameter in the 
ODB, a callback routine in the slow control front-end is 
called, which in turn writes the new value to the associated 
slow control device. By dumping the ODB into a snapshot 
file, the current status of the whole experiment can be cap- 
tured at a specific time, and be reconstructed later from that 
dump. 

In addition, the slow control front-end generates events 
from the slow control parameters and injects them periodi- 
cally via the buffer manager into the event stream written to 
disk for later analysis. 

70.2.2 MIDAS Custom pages 

The whole mhttpd web interface is non-specific, meaning 
that the same pages are used for different experiments. This 
is useful for monitoring the general status, starting and stop- 
ping runs and changing ODB variables. There is however the 
need for experiment-specific pages such as the one shown in 
Fig. [86l Such pages condense information from one or sev- 
eral subsystems into a single display and provide the shift 
operator with a quick overview of a part of the experiment. 
To realise these so-called "MIDAS custom pages", three tech- 
niques are implemented into the mhttpd web server: 

- The HTML language is extended by an <ODB [path] 
[flags] > tag. When the mhttpd process parses the user- 
written HTML file for the custom page, the ODB tag 
is replaced by the actual value of the parameter stored 
in the ODB indicated by its path. Optional flags allow 
various number formats and the possibility to change the 
ODB value though the web interface. 

- Dynamic GIF images can be created by mhttpd. They 
consist typically of a static background image combined 
with overlays of labels, bar graphs or coloured areas based 
on simple script statements. The image is created in mem- 
ory using the GD library [73] and directly streamed to 
the browser for maximum efliciency. The right side of 
Fig. |86]is based on this technique. 

- An Asynchronous JavaScript and XML (AJAX) library 
is available which contains functions for asynchronous 
communication with mhttpd. These functions can be used 
to change parts of a web page using the Document Ob- 
ject Model (DOM) in the background without the need 
to reload the complete page. Functions are implemented 
for reading and writing ODB values, to start/stop a run 
and to generate alarm and system messages. 



The combination of the three methods has proved to be 
extremely powerful in designing web pages which almost 
look and behave like a dedicated program, but have the ad- 
vantage that only a web browser is required to use them. 



Fig. 86 Custom web page of the LXe detector cryogenic system. The 
plots on the left are taken from the history system. 

The usage of custom pages became a very efficient tool 
to control and monitor the complete experiment with only 
a web browser This was a prerequisite for running remote 
shifts, which were introduced in 2011 in the MEG experi- 
ment. 

11 Simulation 

This section provides an overview of the packages used in 
the Monte Carlo (MC) simulation of the MEG experiment. 

The MC simulation is based on the physics event sim- 
ulation (megeve), the detector simulation (gem), the read- 
out simulation and the event mixing (megbartender). These 
packages make use directly or indirectly (via megdb2cards) 
of the MySQL ||74l based Data Base (DB) (more details are 

in mm)- 

11.1 The simulation architecture 

The simulation architecture was designed with a number of 
constraints from the time schedule, manpower and techni- 
cal skills available at the start of the project. Therefore the 
physics and detector simulation packages were written in 
FORTRAN 77, while megbartender was written in C-i-h. 

The FORTRAN 77 packages were organised following a 
modern programming paradigm, i.e. to use an Object-Oriented 
(OO) approach organised in a custom-designed framework 
called rem. The basic unit manipulated by the framework is 
a Module, corresponding to an OO class. 

Modules can be either sub-detector simulation sections 
or service tools such as graphics or input/output. The megeve 
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and gem packages are steered through cards read by the FF- 
READ package. 



11.2 Physics event simulation 

In a Module, megeve, several event types of relevance for 
the experiment can be generated. A non-exhaustive list in- 
cluding exclusive and inclusive decays is (particles in brack- 
ets ai-e not generated) 
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When appropriate the generation includes first order ra- 
diative correction and /u^ polarisation effects [601 . 

The yU^ decay vertices may be forced to be in the target 
or each fi* is propagated through the beam line starting at the 
last quadrupole of Triplet 11 (see Fig. |2]l 5647 mm upstream 
of the centre of COBRA to the decay vertex. 



11.3 Detector simulation 

The beam propagation, interaction in the target and propaga- 
tion of the decay products in the detector are simulated with 
a FORTRAN 77 MC program (gem) based on the GEANT3 
package |77|, which describes the sub-detector responses up 
to the read-out. 

The program is heavily modularised using the rem frame- 
work. That allows easy addition of new Modules for imple- 
menting the various parts of the simulation. Settings depen- 
dent on the run-period are copied from the MySQL database 
into the FFREAD cards, and used in gem. The output is writ- 
ten in ZEBRA format [781 and the code for dealing with the 
ZEBRA package is automatically generated starting from a 
bank description in DZDOC format. 

Each sub-detectors is described by a Modules: DCH, 
Longitudinal TC, Transverse TC and LXe detector 

Additional Modules, e.g. COBRA, gas volume, fi target, 
upstream and downstream elements outside COBRA, com- 
plete the description of the material distribution. 




Fig. 87 A simulated event /j+ —> e+7. 



The GEANT3 routines simulate the energy release in the 
detectors, while the conversion into electrons/photons and 
their propagation up to the read-out devices are managed in 
gem user routines. 

The LXe detector Module includes a very detailed sim- 
ulation of the propagation of each scintillation photon in the 
detector volume, accounting for reflection, refraction and 
polarisation effects. The output is the photoelectron current 
on the PMT photocathodes. In the TC Modules (longitudinal 
and transverse sections) the photoelectron currents on PMT 
(APD) photocathodes are calculated with an analytical ex- 
pression for propagation of the scintillation photons. In the 
DCH Module the GEANT3 energy hits are split into clus- 
ters along the track. Then the electrons are drifted toward 
the sense wires following the position-time relation obtained 
with GARFIELD |16|, including diffusion effects. The out- 
put is wire currents. Figure. [87] displays a simulated signal 
event. 

The gem package stops short of simulating the read-out, 
which is performed by Bartender. 
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11.4 Bartender 



Between the detector simulation and the reconstruction pro- 
gram an intermediate program, called Bartender, is required 
for the processing of MC data. 

The Bartender serves different roles: 

- Conversion of ZEBRA files into ROOT |79 | files 

- Readout simulation 

- Event mixing 

11.4.1 Format conversion 

Bartender reads several ZEBRA files written by gem and 
converts them into ROOT format. Simulation-specific data 
such as kinematics of generated particles, true hit informa- 
tion etc. can be streamed to a "sim" Tree in separate ROOT 
files for further studies. 

11.4.2 Readout simulation 

Bartender simulates the detector read-out electronics and 
makes waveforms. For example, for the LXe detector, wave- 
forms are made by convolution of the single photoelectron 
response of a PMT with hit-time information of scintillation 
photons simulated by gem. PMT amplification, signal atten- 
uation, saturation of the read-out electronics, noise etc. are 
taken into account. Simulated waveforms are encoded in the 
same manner as in the experiment and written to a "raw" 
Tree in ROOT files. 

11.4.3 Event Mixing 

Bartender makes a mixture of several sub-events; rates of 
each event type are set with a configuration file. To study 
the combinatorial background events, sub-events are mixed 
with adjustable relative timing with respect to each other and 
with respect to the trigger; random and fixed timing can be 
selected. This allows simulation of many different pile-up 
configurations with a limited number of event samples sim- 
ulated in the detector 



11.5 DataBase 

Run-dependent parameters such as geometry, calibration and 
analysis constants are stored in a relational database based 
on MySQL (DB). The DB is used by the DAQ front-end, 
analysis and simulation. For simulation, the dedicated pro- 
gram megdb2cards reads the DB and writes the FFREAD 
cards required by gem. Therefore the DB is used consis- 
tently by all packages. 
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Fig. 88 Nonnalised momentum spectram measured in the spectrome- 
ter for data and MC simulation. 
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11.6 Comparison with data 

Figure[88]shows the normalised momentum spectrum of the 
measured in the spectrometer for data and MC events. 
Figure |89(a)||89(b)| shows the hit distribution on the longi- 
tudinal Timing Counter for data and MC events. Figure |90] 
shows the inclusive y-ray spectrum measured in the LXe de- 
tector for data and MC events. The MC events are separated 
into different sources. 
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Fig. 90 Inclusive y-ray spectrum in LXe. Data (dots) and MC (lines). 
MC spectrum is separated into contributions from AIF, RMD and pile- 
up from multiple y-rays. 



Table 9 Resolutions expressed in Gaussian standard deviations (core 
or single side, when required) and efficiencies for all observables: com- 
parison between expected and obtained. 



Variable 


Design 


Monte Carlo 


Obtained 


Resolutions 


Positron(e) 


o-£. (keV) 


200 


315 


306 


o-0„,e„ (mrad) 


5(<^e),5(f^e) 


8(</.e),9(9e) 


9(</'e),9(ee) 


o"z,,v„ (mm) 


1.0 


2.9(z,)/I.0(y,) 


2.4(z,)/I.2C>g 


o"rc (ps) 


50 


65 


102 


Photon (7) 


O-E, (%) 


1.2 


1.2 


1.7 


0"r, (ps) 


43 


69 


67 


(mm) 


4 


5 


5 


o-„.^ (mm) 


5 


6 


5 


Combined(e-7) 


(ps) 


66 


95 


122 


o-Q^^ (mrad) 


11 


16 


17 


Efficiencies 


(%) 


90 


40 


40 


(%) 


60 


63 


63 


e„j (%) 


100 


99 


99 



12 Summary 

The MEG detector described in this paper has been opera- 
tional since 2008. In the period 2009-201 1, it has collected 
~ 3.63 X 10'^ muon decays. 

In the Table |9] the design, Monte Carlo and obtained res- 
olutions are summarised. The obtained resolutions are, for 
some observables, worse than the design values. The differ- 
ences are, at large extent, understood and the Monte Carlo 
resolutions match quite well the obtained ones. 

The influence of the differences in resolutions on the 
analysis sensitivity is approximately evaluated firstly cal- 
culating the number of expected signal events for a given 
branching ratio (B) 

A^,;,, oc ey X fe, (23) 

from Table|9] the decrease in A^.,,,, is ~ 2 due to e^. Secondly, 
the number of accidental coincidences is 

T^acc ^ (r\ X a-E, X cr,^^ X 0-2 . (24) 

This number increases a factor of ~ 10 compared to the 
design value. Those effects combined implies that the sen- 
sitivity of the present MEG experiment is ~ 5 compared to 
the design value. 

The results from the analysis of data collected up to 2010 
are already published and set the most stringent limit to date 



on the existence of the fj.'^ s^y decay. Data collected in 
2011 with statistics comparable to the 2009-2010 sample 
are being analysed. 

The experiment run in 2012 collecting a data sample 
comparable in size to the 201 1 sample. 
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